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Numerous  approximations  have  been  necessary  to  facilitate  evaluation  of 
some  of  the  multiple  integrals;  to  what  degree  the  sufficient  conditions  cited 
above  can  be  relaxed,  without  violating  the  conclusions,  is  unknown. 
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INTRODUCTION 


Detection  of  a  tone  of  known  frequency  and  limited  duration.  In  the 
presence  of  noise,  is  often  accomplished  by  passing  the  received  waveform  through 
a  matched  filter  and  comparing  a  sample  of  the  filter  output  with  a  threshold. 

(For  unknown  phase  of  the  tone,  it  is  the  filter  output  envelope  that  is  compared 
with  a  threshold.)  The  bandwidth  of  this  receiving  filter  is  very  narrow,  in  fact 
approximately  equal  to  the  inverse  of  the  signal  duration;  in  this  manner,  the  noise 
passed  by  the  receiving  filter  is  greatly  attenuated,  and  decent  decisions  about 
signal  presence  or  absence  can  be  made.  Typically,  the  ratio  of  received  signal 
energy  to  noise  power  density  level,  of  the  order  of  10,  is  required  for  low  false 
alarm  probabilities  and  acceptable  detection  probabilities. 

However  if  the  tone  is  phase- or  frequency-modulated  in  a  random  fashion, 
the  spectrum  of  the  received  signal  is  spread  significantly  beyond  the  inverse  signal 
duration,  and  any  receiving  filter  must  be  broadened  to  accept  the  received  signal. 

As  a  result,  considerably  more  noise  is  passed  to  the  receiving  filter  output, 
adversely  affecting  decision  capability.  The  fact  that  the  modulation  is  random, 
that  is,  unknown  to  the  receiver,  prevents  much  coherent  processing  of  the  received 
signal,  if  the  modulation  is  significant.  If  the  deviation  and  bandwidth  of  the 
frequency-modulating  process  were  known  (or  better,  if  the  spectrum  of  the  frequency- 
modulating  process  were  known),  the  question  arises  as  to  just  how  much  could  be 
achieved  in  detection  capability  by  taking  advantage  of  this  information. 

We  will  address  the  problem  of  optimum  detection  of  such  signals  in 
the  following,  with  particular  emphasis  on  the  case  where  the  product  of  observation 
time  (signal  duration)  and  signal  bandwidth  is  much  larger  than  unity.  For  long 
duration  signals,  this  situation  can  easily  arise  with  fairly  small  amounts  of 
phase-  or  frequency-modulation.  We  will  also  consider  simpler  sub-optimum  processors 
and  attempt  to  deduce  their  performance.  The  severe  analytic  problems  preclude 
complete  solution,  and  some  considerable  simulation  effort  will  be  necessary  in 
the  future  in  order  to  give  quantitative  comparisons.  Some  related  work  is  given 
in  Ref.  1  for  the  first-order  term  in  a  series  expansion  of  the  optimum  processor, 
useful  for  low  input  signal-to-noise  ratios. 


LIKELIHOOD  RATIOS 


General  Case 

In  an  observation  of  M  samples,  either  the  condition  noise-alone  or 
signal-plus-noise  prevails.  On  the  basis  of  these  M  samples,  an  optimum  decision 
about  signal  presence  or  absence  is  to  be  made.  The  signal  samples  are  denoted 

by  S  *  (s,  s«),  (1) 

and  the  observation  is  denoted  by 

X  '  (*.  *>'  >0  (2) 

If  we  let  Ho  denote  the  signal -absent  hypothesis,  and  H,  denote  the 
si gnal -present  hypothesis,  the  optimum  decision  rule  is  based  upon  a  comparison 
of  the  likelihood  ratio  (LK)  with  a  threshold  (Refs.  2,  3): 

l*(X)  -  •  (3) 

where  the  noise-alone  probability  density  function  (PPF)  of  observation  X, 

and  Pl(X)  is  the  FDF  of  X  with  signal  present.  We  can  express  (3),  for  any 
signal  statistics,  as 

fdS  p.(X|s)  p,(S) 

~~W)  (4) 

where  p,(X|5)  is  the  conditional  PDF  of  X,  given  signal  values  S,  and  pt(S)  is  the 
a  priori  Pt>F  of  5. 

A  modified  notation  for  (4)  is  adopted  in  the  following,  namely 

L*W-  -  <u«(x)  s)>.  «» 

where  the  ensemble  average,  <>,  is  over  the  signal  statistics,  whatever  they  are, 
considering  observation  X  as  a  fixed  quantity.  (See  also  Ref.  3,  page  132,  eq. 
(1.4)).  (5)  is  an  average  over  the  conditional  LR. 

For  continuous  observation  of  waveform  *(fc)  over  a  time  interval 
where  the  signal  exists  (if  present),  the  IR  in  (5)  becomes 


♦Integrals  without  limits  are  over  the  range  of  non-zero  integrand. 
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u(xR):  = 


(6) 


where  the  ensemble  average,  <  >,  is  over  the  random  signal  process  sfO>  considering 
observation  x(t)  as  a  fixed  waveform.  This  interpretation  of  the  bracket  <>  will 
be  used  throughout  the  following  report. 

Additive  Signal  and  White  Gaussian  Noise 

The  results  above  pertain  to  any  signal  and  noise  model.  We  now 

* 

specialize  to  additive  signal  and  noise.  In  particular,  for  flat  low-pass 
Gaussian  noise  of  bandwidth  K  H?  double-sided),  noise  samples  taken 

every  seconds  apart  are  independent;  then  (for  real  processes) 


■  S\ 


_J _ 


(7) 


and 


(8) 


where  ^ 


is  the  noise  variance,  presumed  known. 


Then  (5)  yields 


uM= 


(9) 


where  the  ensemble  average  is  over  only. 

Since  <T^=N#Wn>  where  is  the  single-sided  noise  density  level 
in  watts/Hz  (presumed  known),  then  in  the  limit  as  the  noise  bandwidth  (9) 

(and  (6))  becomes 


L?  ( xft)  s  10  - 


(10) 


(See  also  Ref.  3,  pages  108-9.) 

If  the  signal  ’It)  and  observation  x(t)  are  narrowband  about  center 


*For  narrowband  signals,  the  assumption  that  the  noise  spectrum  is  flat  in  this 
narrow  band  is  not  very  restrictive. 
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frequency  f#i  and  if  T£»l,  where  observation  time  (signal  duration) 

then  (10)  can  be  expressed  as  (see  also  Ref.  3,  page  160,  eq.  (5.2)) 

Lr(xIH;<+<  *i)  =  OD 

where  Re  denotes  the  real  part,  and  2.10  and  x (0  are  the  complex  envelopes  of 
signal  and  observation  respectively: 

5fc)r  Rt.[slO^f02ir<-t-)]  ,  x/f)  «■  Re  .  (12) 

It  should  be  noted  that  although  we  have  presumed  T(j  »  I,  we  have  not  fixed 
TW5,  the  product  of  observation  time  (signal  duration)  and  signal  bandwidth. 

Also  we  have  not  yet  specified  the  signal  statistics  in  any  way. 

Fixed  Received  Signal  Energy 

A  special  case  of  the  above  is  afforded  when  the  received  signal 
energy  in  the  observation  interval, 

A  & 

B  *  JT  s*/f)  =  ■£■ J  <Hr  |*l»|  ,  (13) 

\  -k 

is  constant,  even  though  the  received  signal  sW  is  a  random  process.  (An 
example  is  the  narrowband  phase-modulated  signal 

sf()=  A  c05[2TT-f,i+  ,  4-a  <  -t<  -tfc  ,  (14) 

where  Tf%»  1>  A  is  a  constant  (but  perhaps  unknown),  0(t)  is  a  random  process, 
and  is  a  random  variable;  the  product  TWS  is  arbitrary.  Then,  E*  /?t]2  ). 

In  this  special  case  of  (13),  (11)  yields 

LK  al0/(tj£>  (15) 

This  simplified  form  for  the  LR  (for  fixed  received  signal  energy)  points  out 

the  magnitude  of  the  analytical  difficulty  of  the  problem  we  are  addressing. 

Namely,  the  ensemble  average  on  the  right-hand  side  of  (15)  is  the  characteristic 
★ 

function  of  the  random  variable 


‘Actually,  it  is  the  moment-generating  function. 
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at  argument  But  (16)  is  a  linearly-filtered  version  of  random 

process  5/f),  with  fixed  waveform  &(t)}  recall  that  aft)  is  considered  fixed 
in  the  ensemble  average  in  (15).  And  since  calculation  of  the  PDF  of  the  output 
of  a  linear  filter  (for  other  than  Gaussian  input  processes)  has  remained  an 
unsolved  problem  in  the  literature  for  many  years  now,  there  is  no  hope  of  exact 
evaluation  of  (15)  for  general  random  processes  Rather,  all  we  can  hope 

for  is  good  approximations  for  some  signal  processes  under  certain  conditions. 

The  situation  for  the  more  general  case,  (11),  is  even  worse,  due  to  the 
additional  randomness  of  the  received  signal  energy.  Notice  that  the  only 
restriction  on  the  signal  statistics  in  (15)  is  that  stated  in  (13);  we  are  not 
restricted  to  the  model  of  example  (14). 

In  the  following,  we  shall  generally  make  the  assumption  about  the 
received  signal  energy  being  constant;  notice  that  this  disallows  sfO  or  i/t) 
from  being  a  Gaussian  random  process.  We  shall  also  concer  rate  on  narrowband 
signals,  for  which  however,  TX  is  arbitrary  initially.  Hence  the 

pertinent  expression  for  the  LR  is  (15). 

It  is  perhaps  worthwhile  at  this  point  to  fix  the  ideas  presented 
thus  far  with  a  couple  of  examples.  In  particular,  averaging  with  respect  to 
signal  statistics,  while  holding  the  received  waveform  (observation)  fixed,  will 
be  illustrated.  The  received  signal  energy  is  not  constant  for  these  two  examples. 
Low  Pass  Example 

Let  the  received  signal  be 

sft)=  A  w(t),  (17) 

where  A  is  random  but  positive,  and  is  a  known  deterministic  real  function. 

Since  the  only  randomness  in  the  received  signal  is  through  the  scaling  parameter  A, 
(10)  becomes  ^  ^ 

U  -  »*] ^ 

where 

\  4 

a*  ^-J*&  xtewlW  ,  dim'd). 


(18) 


(19) 
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Notice  that  the  observation  xft)  appears  in  the  LR  only  through  the  coherent 
integration  yielding  parameter  a.  Now  we  can  express  (18)  as 

Lfi  =  I  dA  j>(A)  e*j>[Aa-Al0.  (2°) 

But  since 

<2-  LK  =  £<M  f(A)  A  e*r[A«  -  A1  k]  (21 ) 

is  positive  for  any  PDF  flA),  when  A  is  limited  to  positive  values,  then  LR 
is  monotonical ly  increasing  with  the  parameter  Q  .  Hence  comparison  of  LR  with 
a  threshold  is  synonymous  with  comparison  of  parameter  a  with  a  threshold. 

Thus  the  optimum  processor  is 

\  H, 

1  *  x(f)«lO  >  AUsUJ,  (22) 

V  H. 

regardless  of  the  PDF  p(X),  when  random  variable  A  is  limited  to  positive  values. 
(22)  is  recognized  as  the  standard  coherent  correlation  detector  operating  on 
observation  x(f). 

Narrowband  Example 

Let  the  received  signal  be 


S  ft)  =  A  fie  ft)  e/xp  (i  2r-fei  +  <  4>)] ,  4  «  *  c  \  >  (23) 


where  A  is  random  and  positive,  and  is  random  and  uniformly  distributed  over 
2-rr.  Variables  A  and  are  independent;  there  is  no  need  to  consider  negative  A, 
since  this  effect  can  be  absorbed  by  values  of  T ■  .e  complex  envelope  is  then 

A  U)  f/0,  (24) 


where  fft')  is  a  known  deterministic  low-frequency  complex  function.  Since  the 
randomness  of  ?l()  arises  only  through  A  and  <{>,  (11)  yields 


=  <^I^Aa)  e*p(- A*l>)^  =  (25) 


6 


where 


A  *  ^  • 

Notice  that  observation  &ft)  appears  in  the  LR  only  through  the  coherent 
integration  yielding  parameter  a.  But  since  (using  (25)) 

^  u  =  £dA  fik)  A  T,(a.)  «f(-  K  i.) 

is  positive  for  any  PDF  pt4),  then  LR  is  monotonically  increasing  with  the 
parameter  a.  Hence  comparison  of  LR  with  a  threshold  is  synonymous  with 
comparison  of  a  with  a  threshold.  Thus  the  optimum  processor  is 


regardless  of  the  PDF  p(A).  (28)  is  recognized  as  a  threshold  comparison  of 

the  envelope  of  the  output  of  a  matched  filter  operating  on  observation  $fl). 


(26) 


(27) 


(28) 
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SERIES  EXPANSION  FOR  LOGARITHM  OF  LIKELIHOOD  RATIO 

The  starting  point  for  this  development  is  (15),  which  applies  to 
arbitrary  narrowband  signal  statistics,  subject  to  fixed  received  signal  energy. 
When  we  define  random  variable 

V  =  {  dt  J  It)  /W,  (29) 

(using  the  fact  that  received  signal  ift)  and  observation  xfi)  are  non-zero 
only  in  tt», "{},)),  then  (15)  can  be  expressed  as 

LK  =  <C^W)>>  (30) 

where  the  ensemble  average  is  over  the  signal  statistics.  Now  from  (29),  we 
find 

<v>  -  M-  Re  Jdi  sit)  <  l*k£>  ,  (31) 

where  <f will  be  called  the  deterministic  signal  component.  Then  define  the 
residual  of  random  variable  v  as 

r  =  v  -  <v>  »  Re  |  ^  sH)  2*lt),  (32) 

where  we  have  also  defined  the  "ac"  component  of  the  signal  complex  envelope  as 

flit)-  sit)-  (33) 

(Later  we  shall  generally  consider  signal  models  for  which  the  deterministic 
signal  component  is  absent,  i.e., 

<  £lt»  -  0,  (34) 

which  is  the  usual  case  for  phase-random  narrowband  waveforms.  For  example,  (24) 
was  just  such  a  case. ) 

Substitution  of  (32)  in  (30)  yields 

=  !r  +  <v^)  . 


(35) 


At  this  point,  we  define  the  characteristic  function  of  random  variable  y 
in  (32)  to  be  the  ensemble  average 


£('?)= 


Once  again,  we  note,  with  reference  to  the  definition  of  r  in  (32),  that  the 
randomness  of  r  is  to  be  considered  due  to  that  of  ac  signal  cj(t)  (through  (33)), 
and  that  observation  must  be  considered  fixed  in  (32)  for  this  ensemble 
average  in  (36). 

The  natural  logarithm  of  the  characteristic  function  of  Y  can  be 
expanded  as  . 


•  i  Or)  ■  £4*. 


where  is  the  k-th  cumulant  of  random  variable  r  (See,  for  example.  Ref. 

4,  Chapter  3).  The  sum  starts  at  k  =  2  since  (see  (32)) 

<r>  =  <V-<v>>  =  0.  (38) 

For  example,  we  have  (Ref.  4,  page  71,  eq.  3.43) 

■X'rt=  <r-> 

7$  -  <rJ> 

X!,'’  --  <r4>  -  3<K>  (39) 

Since  the  natural  logarithm  is  a  monotonic  transformation,  comparison 
of  LR  in  (35)  with  a  threshold  is  synonymous  with  comparison  of  In  LR  with  a 
(different)  theshold.  Then  using  (35),  (36),  and  (37),  we  find 

In  LR  =  - Kv>  +  In  tr(i) 

-  +  <*>  _  it  *  Jr'*" (40) 

where  is  the  k-th  cumulant  of  random  variable  v;  see  (32).  The  first  term 

in  In  LR,  namely  -E/hti  is  a  constant,  independent  of  the  observation  i(f)>  and 
can  therefore  be  ignored  (i.e.,  absorbed  in  a  modified  threshold  in  the  decision 
comparison). 

The  second  term,  <v>,  in  |w  LR  is  proportional  to  a  coherent 
correlation  of  the  received  waveform  with  the  deterministic  component  of  the  signal 
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waveform;  consider  the  following  development  (where  co*  *  2n-f#)  ; 

xW<»W>  =  jj-JiJt  K«(x»='At}Oe^(t)e"4+]> 

=  -0*  feJjBe'*4] 

=  "j^  Re  f  A  it)  <'s*lt))>  »<v>,  (41 

where  we  have  used  (12),  the  fact  that  t4j»|>  and  (31).  (As  noted  in  (34), 
this  term  will  generally  be  zero  for  our  signal  models  of  interest;  that  is,  there 
will  be  no  deterministic  signal  component  in  the  received  waveform,  and  hence  the 
<v>  =  term  will  be  absent  from  the  In  LR  expansion  in  (40).) 

All  the  other  terms  in  In  LR  in  (40)  involve  r  defined  in  (32)  (or 
v  defined  in  (29)).  In  particular,  since  (32)  can  be  written  as 

r  -  (42 

then  we  have,  using  (39),  the  second  cumulant 

y?  =  02>  -  ^  If* 

-  A,  &K) IL) ^ ? If,)  9* M.)^>  +  ComjA*  Co*  (43 

where  we  have  used  the  fact  that 

<  g  It)  t  tt.)>  -  0 ,  o)U,t  (44 

for  complex  envelopes  of  narrowband  random  processes;  see  for  example.  Ref.  3, 
page  53,  eq.  (5.12).  Now  it  may  easily  be  shown  that  each  of  the  two  terms  in 
(43)  is  real;  therefore 

T^-  ik  &  sfc)  R* M,)>  (4s: 

where  is  the  autocorrelation  of  ac  component  off),  assumed  stationary  in  the 
observation  interval;  that  is, 

<  afe)  s%)>  =  **  k-  i,  k  £  fe,  h).  (46: 


An  alternative  expression  to  (45)  will  be  developed  below. 


The  remaining  cumulants  in  (40)  require  higher-order  moments 
of  aft),  and  must  await  evaluation  until  we  specify  the  exact  signal  model  of 
interest,  which  we  have  not  yet  done.  However  it  is  very  important  to  observe 
that  if 

(47) 

were  a  Gaussian  random  variable,  then  all  its  cumulants  ,  for  would 

be  identically  zero,  and  the  series  for  In  LR  in  (40)  would  terminate  with  the 
term  (45).  Since  xK)  in  (47)  is  to  be  considered  fixed  in  this  particular 
statistical  consideration  of  r,  it  is  seen  from  (47)  and  (33)  that  r  is  the  sum 
of  a  large  number  of  statistically  independent  components  if  the  time-bandwidth 
product  of  the  received  signal  is  large,  that  is,  if 

TWj  I,  (48) 

and  if  alt)  is  fully  random  (to  be  explained).  This  is  due  to  the  fact  that  the 
integral  in  (47)  is  over  a  time  duration  of  T  seconds  and  that  &K)  has  a  statistically- 
independent  "wiggle"  approximately  every  l/Ws  seconds.  Then  by  appeal  to  the  Central 
Limit  Theorem,  r  will  be  nearly  Gaussian  when  (48)  obtains,  and  the  cumulants  for 
k*3  in  (40)  can  then  be  well  approximated  by  zero.  We  will  attempt  to  quantitatively 
justify  this  claim  later,  for  a  particular  signal  model. 

The  expression  for  the  1m  LR  is,  from  (40), 

In  LR  =  -ft  +<v>  +  +  (49) 

•vW 

Now  /Ch  is  actually  a  random  variable,  governed  by  the  statistics  of  the 
observation  X  K)  >  "t*  ;  remember  we  have  already  averaged  over  the  statistics 

of  s  It),  in  (36)  and  (37).  Furthermore  has  different  statistics  for  signal- 

present  (H,)  versus  signal-absent  (H»).  Although  the  sum  of  the  random  variables 
in  (49)  is  to  be  compared  with  a  threshold,  it  is  not  the  absolute  level  of  the 
mean  values  of  ^/k!  >  under  H»  and  >  that  is  important;  rather,  it  is  the 
di fference  in  these  means  under  H,  and  H,  that  is  important.  After  all,  means 
can  be  absorbed  in  a  modified  threshold.  Therefore  we  will  attempt  to  determine 
the  means  and  standard  deviations  of  some  of  the  random  variables 
and  make  some  quantitative  statements  as  to  when  and  where  (49)  can  be  terminated. 


Another  important  point  should  be  made  here  regarding  (49).  One 
sub-optimum  processor  is  afforded  by  simply  terminating  (49)  with  the  term  ytC/2- 
Now  even  though  the  rest  of  the  terms  in  the  series  may  not  be  small,  and  therefore 
cannot  be  dropped  from  the  optimum  processor,  it  is  conceivable  that  the  decisions 
yielded  by  this  sub-optimum  processor  often  agree  with  those  yielded  by  the 
optimum  processor  (each  processor  with  its  own  threshold).  This  will  be  the  case 
when  (48)  is  true;  an  important  question  is:  how  small  can  TV^  get  before  the 
optimum  processor  signi ficantly  outperforms  this  sub-optimum  processor?  An  example 
in  appendix  A  shows  no  difference  in  performance  for  very  small  TWS<>  namely 

T  <  »/Ws  .  Thus  at  the  two  extremes  of  very  large  and  small  this 
sub-optimum  processor  performs  just  as  well  as  the  optimum  processor. 

It  is  worthwhile  to  note  that  if  we  were  to  expand  the  LR  in  (35) 
in  a  series,  we  would  obtain  a  far  less  useful  result;  namely 

LK  =  SSKf  4-  <»)  ^  f‘k)>  (50) 

Now  this  series  would  not  terminate,  even  for  r  equal  to  a  Gaussian  random  variable. 
And  for  r  near-Gaussian,  it  is  not  obvious  how  many  terms  of  (50)  need  to  be 
retained  for  a  good  approximation.  Thus  (49)  is  expected  to  be  more  useful, 
especially  when  (48)  is  true. 

Example  of  Exact  LR  for  Several  Narrowband  Pulses 

By  way  of  illustration  of  some  of  the  techniques  above,  consider  the 

signal  model 

s  ft)  -  'Ke  \  5  ^  *?['  *4*  +  •  4^  ;  ^  *  h  >  (51 ) 

where  are  independent  random  variables  uniformly  distributed  over  2rr,  l&lt)] 

are  known  complex  deterministic  low-frequency  waveforms,  are  known  amplitudes, 

and  are  known  center  frequencies.  We  assume  that  the  complex  waveforms 

^5* (B  «Ap 0 are  orthogonal,  perhaps  through  time  or  frequency  separation, 
and  that  T^»l,  where  T*  is  the  duration  of  5  •  Then  it  follows  that 

the  received  signal  energy,  J £  is  virtually  independent  of  {fa] ,  and 

hence  is  non-random. 

The  signal  complex  envelope  is 

sit)  -  +  ,  (52) 
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for  which,  obviously,  the  deterministic  component  “0.  Then  (32)  yields 

*  s  nT  Kt  \ 2$ It)  ^4 5 *lb)  e*r [->' k  - i  **)*] 

The  ensemble  average  in  (35)  is  then 

<(^ry  =  TT 2$  «*/>[-'$]] 

-  TT  r0  (j^  |  f  *  ft)  f/fd i ( ) 


Hence  (35)  yields  (with  the  aid  of  (31)) 

U  LR  =  *  jf  ^  )  I#  (55) 

A  block  diagram  of  the  optimum  processor  for  this  signal  model  is  shown  in 
figure  1  below.  The  non-linearity  In  T^C  )  is  approximately  a  squarer  for  small 
arguments  (inputs),  and  is  approximately  linear  for  large  arguments.  Thus  coherent 
processing  of  each  component  is  accomplished,  and  then  followed  by  non-linear 
envelope  detection  and  summation. 


Fit 


iwre 


Op  Vi  wuw 


Processor  -|&r  Several  Narrou>kai»4  Pulses 
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BLOCK  DIAGRAM  FOR  PROCESSOR  YIELDING 


In  this  section  we  will  give  a  block  diagram  of  a  processor  which 
accepts  complex  envelope  waveform  ■£*<£<■£»>  at  its  input,  and  emits 

the  quantity  y^(2  To  do  this,  we  first  define  correlation 

K Ir)  =  fr) -  R* (*)  ;  b)=o.  (56) 

An  example  of  a  signal  process  for  which  K (even  though  the  deterministic 
signal  component  <3r(t)>  =  o)  is  given  in  appendix  B.  Furthermore,  it  is  shown 
there  that  the  deterministic  signal  component  corresponds  to  pure  tones 
in  sft)  of  known  phase,  whereas  the  ac  component  with  corresponds  to 

pure  tones  in  s(t)  of  random  phase.  The  quantity  must  be  real,  since 

it  gives  rise  to  an  impulse  in  real  power  spectrum  at  f=o. 

Substitution  of  (56)  in  (45)  yields 

ilC  ‘  'Si, 4 xfcl&RM,).  (57) 

This  processor  is  indicated  in  block  diagram  form  in  figure  2.  The  power 


2.  Processor-  Generate 
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spectrum*  £0K)  is  defined  as  the  usual  Fourier  transform  of  correlation 


(-f)  =  Jdr  (~>  f r)  \  k),  (58) 

and  H(f)  is  a  low-pass  Eckart  filter. 

The  way  to  show  that  the  lower  branch  of  figure  2  generates  the 
second  term  in  (57)  is  as  follows:  the  output  of  the  bottom  integrator  in 
figure  2  is 

|  f  4w h^'u)a(«)|  =  {{in  ilu)s*(v)  J dt  h  w) .  (59) 


But 

=  (  #  exj>[i  2^rf  Cv-u)] 

=  f #  ' 

where  h  and  H  are  a  Fourier  transform  pair  (impulse  response  and  transfer 
function),  and  we  used  (58).  Thus  the  last  quantity  in  (59)  becomes  precisely 
the  last  term  in  (57). 


Since  fr0lf)  =  W  must  be  non-negative,  since  (ct  (f ) 

is  a  legal  spectrum.  Therefore  the  choice  of  )H(f)|2  in  figure  2  is  always 
legal . 
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The  upper  branch  of  the  processor  in  figure  2  depends  on  the 
statistic  \s*m'  (see  (57)),  which  corresponds  to  a  coherent  cross¬ 
correlation,  over  the  entire  observation  interval,  of  the  received  waveform  with 
a  constant  carrier  at  i ,,  followed  by  square-law  envelope  detection  at  the  end 
of  the  observation  interval.  Whether  this  coherent  component  is  significant 
depends  on  the  relative  sizes  of  flgfc*)  and  see  appendix  B  in  this 

regard.  It  should  be  made  clear  that  the  upper  branch  of  figure  2  is  present 
even  when  the  received  signal  has  no  deterministic  component,  i.e.,  when 
<(sK:X>-0.  The  distinction  between  the  two  cases  is  well  illustrated  by  the 
special  case  of  a  pure  tone  signal  of  random  phase,  cos This 
signal  has  no  deterministic  component  because  the  average  over  gives  zero. 
However,  cross-correlation  of  this  signal  with  a  local  carrier  at  -f,,,  with 
in-phase  and  quadrature  components,  followed  by  squaring  and  summing  outputs, 
affordj a  significant  amount  of  coherent  signal  processing  gain  that  cannot  be 
disregarded;  in  fact,  the  performance  is  but  slightly  poorer  than  detection  of 
a  signal  with  known  phase,  4?  (see  Ref.  3,  pages  88  and  155).  When  the  signal 
satisfies  both  of  the  properties 

S  =  0  QiJ  (*)  "  0,  (61 ) 

we  then  say  that  the  signal  is  fully-random;  this  case  will  occupy  most  of  our 
attention  in  the  following. 

The  lower  branch  of  the  processor  in  figure  2  is  a  linear  Eckart 
filter  followed  by  an  energy  detector.  Since  the  linear  filter  impulse  response 
is  non-zero  in  duration  (being  the  Fourier  transform  of  H(Q),  the  finite- 
duration  input  2ilO  is  smeared  by  this  filtering  operation;  hence  the  integration 
in  this  lower  branch  is  over  all  time,  not  just  as  in  the  upper  branch. 

This  Eckart  filter  followed  by  an  energy  detector  is  exactly  the  filter 
characteristic  which  maximizes  the  output  deflection;  see  appendix  C. 

Furthermore,  the  linear  filtering  action  of  H(f)  in  figure  2  is 
such  as  to  perform  coherent  integration  on  the  random  input  signal  over  as  long 
a  time  as  possible.  To  see  this,  we  notice  that  the  output  of  filter  H(f)  can 
be  written  as 

-  j<hr  ilr)  Kdr-r).  (62) 
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But  since  the  passband  of  H(f)  Is  the  same  as  the  bandwidth  of  (r,{F),  the 
impulse  response  duration  of  Mt)  is  approximately  the  correlation  time  of 
■R0h).  That  is. 


(2*r 


(63) 


transforms  into 


Vilu)^(u-r) 


KM 

W 


(64) 


Thus  the  effective  integration  time  in  (62)  is  approximately  the  correlation  time 
of  the  random  signal  process,  and  is  as  long  a  time  as  coherent  integration  could 
reasonably  be  anticipated.  The  squarer  and  integrator  following  H(f)  in  figure 
2  can  be  looked  upon  as  a  sum  of  squared  envelopes,  for  the  continuum  of  possible 
relative  time  delays  that  the  filter  impulse  response  and  input  waveform  can 
take  on  (see  (62)).  For  a  very  narrowband  spectrum  the  correlation  time 

of  will  be  large,  and  considerable  coherent  processing  is  achievable  via 

(62),  prior  to  energy  detection.  This  processing  of  the  random  signal  is 
achievable  even  if  the  upper  branch  of  figure  2  is  absent;  that  is,  even  if  there 
are  no  pure  tones  of  random  phase  in  the  signal. 

An  alternative  interpretation  of  the  quantity  in  (67)  is 

possible.  Let  the  voltage  density  spectrum  of  the  received  waveform  complex 
envelope  be  defined  as 

Xlf)=  e/xp  (65) 


This  is  fine-grained  frequency  analysis,  i.e.,  resolution  >/r  in  frequency.  Then 
it  is  easily  shown  that  (57)  can  be  written  as 


K)' 


IxwT+i# 


(66) 


Realization  of  this  form  requires  a  fine-grained  frequency  analysis  over  iio,W, 
weighting  of  the  energy-density  function  |X(f)|*  with  an  Eckart  filter  function, 
and  summing;  however  the  weighting  at  zero- frequency  is  distinctly  different 
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from  that  at  non-zero  frequency.  Namely  the  weighting  is  (see  (45)  and  (56)) 


(67) 


Since  the  expression  for  7%  in  (45)  involves  Hal tj-i,) 
only  for  the  range  of  arguments  i*  <  tr,,  -ta  <  -fc>  ,  it  is  seen  that  only  the 
values  of  RjJtj)  for  |u|<  T  are  relevant.  The  question  then  arises  as  to 
whether  the  procedure  in  (56)  et  seq.  is  necessary.  This  problem  is  addressed 
in  appendix  D,  where  it  is  shown  that  other  approaches  may  not  yield  legal 
choices  for  |H|f)|a;  that  is,  |H(f)|*  may  be  required  to  be  negative  for 
some  values  of  -f,  which  is  impossible.  At  least,  (56)  guarantees  a  legal  choice 
for 
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MEAN  VALUES  OF  X 


It  was  noted  in  (49)  et  seq.  that  cumulants  iXj?}  are 
actually  random  variables,  since  they  are  functionals  of  the  observation 
In  an  attempt  to  determine  what  terms  in  (49)  are  most  important,  we  begin  by 
first  evaluating  the  mean  values  of  X  >  both  with  and  without  signal  present. 
For  hypothesis  (noise-only),  the  mean  value  of  is  available  from 

(45)  as 

_  ^  _ 

-L  x%)  =  ~r  JJ&,  ^  \  M-,), 

1, 


(68) 


*  . 

where  an  overbar  denotes  an  ensemble  average  with  respect  to  the  statistics 
of  the  received  waveform.  For  the  complex  envelope  noise,  we  have  (Ref.  3, 
ch.  2,  (4.10)  and  (5.15)) 


Substitution  in  (68)  yields 

_ _  A  A 

-  iir  j  * r.(o)  -  ^7 1  *  <iiW> 

<£*U»r>. 


Now  if  < I tt)>  =  0,  no  deterministic  signal  component,  which  is 
the  usual  case  for  phase-random  narrowband  signals,  and  the  only  case  that  we 
will  consider  henceforth,  then  g  lfc)  =  Iff),  and  we  have,  using  (13), 

*  %  <£*uer>  *  i  <»>  =  i  ■  (7 

*This  average  is  over  an  ensemble  of  received  signal  and  noise  realizations  which 
are  both  totally  independent  of  the  earlier  signal  ensemble  average,  although  the 
same  signal  statistics  are  valid. 
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(From  (70)  and  (71),  we  extract  the  relation  (to  be  used  often  later) 

TU.Io)  -  rSs  lo)  =  IB.)  (72) 

Since  the  quantity  e/aJ,  must  be  large  compared  to  1  in  order  to 
get  decent  detection  performance,  even  for  coherent  processing  of  a  known 
deterministic  signal,  it  is  seen  that,  in  this  more-random  signal  case,  the 
mean  value  \  ^fU))  is  very  large,  even  when  no  signal  is  present.  However, 
of  more  importance  is  the  di fference  of  means  of  ^  ,  with  and  without 

si gnal . 

For  signal  present,  (45)  yields  (since  signal  and  noise  are 

independent) 

iXW)  =  -*T  $  +  sUrt)]  \  fe-C).  (73) 

"k a 

Now  recalling  =  we  have 

sit,)  =  <  S  It,)  S*fi£>  =  K  +>)  =  K  Q,  (74) 

and  (73)  becomes,  upon  use  of  (69)  and  (72), 


where  we  have  defined 


* 

See  footnote  on  previous  page. 


20 


The  difference  of  mean  outputs,  with  and  without  signal  present, 
is  available  from  (71)  and  (75): 


The  size  of  this  difference  depends  critically  on  the  value  ofVThe  quantity 
li'  is  a  measure  of  the  effective  number  of  statistically  independent  wiggles 
of  sH)  in  see  Ref.  5,  pages  5-10.  Its  value  is  strongly 

dependent  on  the  size  of  Ra (»).  For  example,  if  Rafc)  does  not  decay  with  T, 
but  remains  constant  at  its  origin  value,  then  ^>lt)«|  and  X*.  =  I ;  this 

is  consistent  with  the  physical  observation  that  an  unmodulated  signal  has 
only  one  statistically  independent  value.  On  the  other  hand,  if  (61)  is  true, 
then  for  T>  3/Ws ,  (76)  yields 


_L  Jeff  If)  _  J_ 

T  '  T ’ 


(79) 


where  we  have  employed  (77),  (74),  Parseval's  Theorem,  and  defined  the  statistical 
bandwidth  of  fully-random  received  signal  as 


(80) 


In  this  case,  Xx<  j  )  generally  X2  can  be  much  smaller  than  unity  if  TV4»I- 
Notice  that  TWS  only  has  to  be  moderately  larger  than  1  in  order  for  (79)  to 
be  a  rather  good  approximation. 

In  general,  (76)  ai.d  (77)  show  that  <  I.  The  nature  of  the 
dependence  of  T*.  on  Rt|«>)  and  T%  can  perhaps  best  be  illustrated  by  an 


Since  Gra(f)  is  the  power  density  spectrum  of  complex  envelope  5^  it  is  a 
low-pass  spectrum  centered  about  f  =  o.  Thus  measures  the  "width"  of  65(f) 
on  both  the  positive  and  negative  frequency  scales  together.  Alternatively,  W* 
measures  the  width  of  the  narrewband  spectrum  of  sfi)  about  its  carrier  frequency. 
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example.  Let  the  signal  correlation  be  exponential: 

(r)  =  [j?,,  (o) -  ?9  h)j  e*p(-  %  M)  +■  \  f*6)  '  (81 ) 

equivalently,  using  (77), 


f 


4-  C 


=  Kh)/Bt  (o). 


(r)  *  (\~  c) 

Then  from  (76), 

Tx*  c2  7c-  >ce' "-(l-c)^],  ^TW„ 


(82) 


(83) 


As  particular  cases. 


li*M  1,.  »  I 

)Ho 


X-**e 


x.-'-jrO-i) 

rx=  I  far  a)| 


-f>V  fill  C  , 

>f  c  *  0  ; 

if 

if  c  *  i  ; 

>f  c  ~  o. 


t\5 


-r  -  jjf-l  +•  e 
^  2X1 


-lx 


C  -  0; 


(84) 


These  relations  agree  with  those  stated  above.  A  plot  of  (83)  is  given  in 
figure  3.  The  smallest  value  of  T»,  and  hence  the  smallest  value  of  the 
difference  of  means,  (78),  is  attained  when  e  =  0,  i.e.,  (-»)  =  0.  .  This 

situation  will  yield  the  most  difficulty  in  signal  detection,  and  is  the  one 
we  concentrate  on  henceforth;  that  is,  we  consider  only  the  fully-random  signal, 
characterized  by  (61).  Notice  from  figure  3,  that  for  c =  0,  X2  is  very  well 
approximated  by  ^|  -  f  for  yWs  as  small  as  1.5.  In  fact,  (TWS)  ' 

is  a  good  approximation  to  I,  for  TX>  10. 

When  the  signal  is  obtained  via  Gaussian  frequency-modulation  of 
a  carrier,  the  pertinent  requirements  on  the  spectrum  6rm(f)  of  the  modulating 
process,  to  yield  a  fully-random  signal,  are  derived  in  appendix  E  and  are 
related  to  the  behavior  of  near  -f  ~0.  This  case  of  Gaussian  frequency- 

modulation  with  property  (61)  will  occupy  the  remainder  of  our  attention. 
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VARIANCES  OF  ^ 

For  noise-alone,  the  mean  value  of  )£?  is  available  in  (71). 

To  evaluate  the  variance,  we  start  with  (45)  and  (74)  and  form  the  mean  square 
quantity 

'  MfJ <*, WV  M.)-  (85) 

Now  since  n (-t)  is  a  complex-envelope  Gaussian  process,  we  can  employ  Ref.  6 
and  (69),  and  obtain  for  the  fourth-order  average  in  (85),  the  value 

*4)  +  &  fe"  .  (86) 

Substitution  of  (86)  into  (85)  yields 

m  s[  ^(4  =  ^  *  W 1°) +■ lii.it-  ^ ].  '87> 

Employment  of  (72)  and  (75)  yields 

=(00+  <88> 


Combining  this  result  with  (71),  the  standard  deviation  of  interest  is 

-  f,/?. 


(89) 


Notice  that  this  quantity  can  be  significantly  smaller  than  (71),  the  mean  value, 
if  T*«  ).  This  latter  situation  will  occur  if  TW* »  I  ;  see  (79). 

We  are  now  able  to  form  a  deflection  statistic  for  the  random 
quantity  in  the  |n  LR.  Using  (78)  and  (89),  it  is  defined  as 


.  E  r-7 

$D *  ** 


(90) 
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This  relation  has  been  derived  without  any  specific  assumptions  about  the 
detailed  signal  statistics,  except  for  the  randomness  assumption  (61).  If 
our  processor  is  limited  to  just  the  jrtf?  term  of  the  In  LR  processor  in 
(40)  (recall  that  the  <v>  term  is  zero  by  (61)),  we  would  require  that  (90) 
be  somewhat  larger  than  unity,  in  order  that  good  performance  be  attainable. 

A  reasonable  ballpark  figure  for  low  false  alarm  probability  (10  )  and  decent 

detection  probability  (.5)  is  to  have  deflection 

4  -  \f>0  i  ( 

higher-quality  performance  will  require  still  larger  values  of  Then  (90) 
indicates  that  we  will  have  to  require 

x  -  jL.  a/ZST  ( 

But  if  TYJj  »  I,  fr/H  is  going  to  have  to  be  rather  large  to  get  good 
performance  from  jX**  alone,  for  this  fully-random  signal.  The  factor  xj4 
is  a  measure  of  the  penalty  of  having  to  detect  a  fully-random  signal  rather 
than  a  completely- known  deterministic  signal,  by  employing  just  the  second 
cumulant,j)£;,  of  the  infinite  series  (49). 

For  the  signal -plus-noise  hypothesis,  the  mean  value  of 
is  available  in  (75).  To  evaluate  the  variance,  we  start  with  (45)  and  (74) 
and  form  the  mean  square  quantity 

**a 

where  fourth-order  average 

FltA.vO 2  [*  M  *  *  U][A)  +  fvj  [1  fy*  4  [A> .  ( 


Using  the  statistical  independence  of  signal  and  noise,  their  zero  means,  and 
property  (44),  we  find  (using  an  obvious  shorthand  notation) 


I 
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(95) 


+  sfs,  1,  <  +  3,  3*  +  w,  »*  w,  n* 


All  the  terms  in  (95)  have  been  encountered  previously,  with 
the  exception  of  the  leading  term.  In  order  to  evaluate  it,  we  will  need  to 
be  specific  about  the  signal  statistics.  As  noted  below  (84),  our  interest 
is  in  Gaussian  frequency-modulation  of  a  carrier,  with  property  (61).  The 
pertinent  second-order  characteristics  of  this  signal  model  are  presented  in 
appendix  E.  In  order  to  evaluate  the  fourth-order  moment  required  in  (95), 
we  employ  (E-6)  in  the  form 


ift)  =  A 


exp 


*»i  (u)  4-  i 


-k  *  ^  K  ^  > 


(96) 


where  r* It?  is  the  Gaussian  frequency-modulating  process,  to  obtain 


c  5*  s  5* 

°i  •> 


4  P  "t 

A  ■&*}>[_(  clu  h(v) 


+  i  2-» 


-  h  ^^\_  \  lir  j"  Ju  m  lu)w(u;  i)] 


(97) 


where 

=  (98) 

and  TJOd  is  the  unit  step  function: 


Equations  (97)  and  (98)  hold  for  any  values  of  4a,  t*  in  ft#,  O* 
irrespective  of  their  location  and  ordering  in  size. 
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Now  since  frequency-modulating  process  m(u)  in  (97)  is  zero-mean 
Gaussian,  the  integral  on  u,  denoted  by  is  also  zero-mean  Gaussian,  and  we 
can  immediately  write  (97)  as 

S',  5?  s,  5  *  =  A  exp [i  2^"]  - 

=  A4  e*p[-2-n-x  Jfyn  Jv  R*  [k-v)  w  (u;  t)  *(v ;  £)) 

--  A+  **p[- 2tt5  J  d-f  £r*  If)  | W(f ;  t)| 2 


=  J  dw  e*f('i2rrfu)w(u;i) 


Here  we  have  employed  the  Fourier  transform  relation  between  frequency-modulating 
spectrum  Gc„(f)  and  its  correlation  R*  h}  (see  appendix  E),  and  used  (98). 

Equation  (100)  is  a  compact  expression  for  the  fourth-order  moment, 
and  can  be  numerically  evaluated  by  one  integral  when  values  for  -(i; 
are  specified.  However,  it  is  not  directly  suitable  for  our  use  in  (95)  and 
(93).  Rather,  we  expand  j'W'P  by  means  of  (101),  and  obtain 

^  b  S3  S*  = 

-  tr* ff, - h)  -  <r7 l-k-  +J  -  <r2 /fr4*)  -  c rVV  ^jjj,  (it 


(ic 


where 


(T1  |r)  s  2.]"  c}f  ^  [l-  Cos(>irfr)] 


In  appendix  E,  this  real  even  function  is  interpreted  as  the  variance 

of  the  "random  walk"  process  2-rrJtlui«ilH)  during  an  Interval  of  |rl  seconds. 
Thus  we  have  a  closed  form  expression  for  the  fourth-order  moment  in  (102), 
provided  that  (103)  can  be  evaluated  in  closed  form. 

Now  also,  for  the  frequency-modulated  signal  process  in  (96),  we 

have 

(-t,-  tj  -  2  ft,)  5^x)  =  A  &KJ>  [l  2-rr  J ^  luj] 


(104) 


which  is  real.  The  last  step  follows  from  appendix  E,  or  by  setting  * 't*. 
in  (97)  and  (102),  and  noting  from  (103)  that  O'3  (o)  *  0.  But  now  (104)  enables 
us  to  write  (102)  as 


Kk-h)K,k,-h) 


(105) 


in  terms  of  second-order  statistics.  This  equation  is  exact  for  signal  model 
(95).  For  comparison,  if  5(0  were  a  Gaussian  process,  the  right-hand  side 
of  (105)  would  have  been  expressible  in  terms  of  second-order  statistics  (according 
to  Ref.  6)  as 


^5  +  h) 


(106) 


(The  approximation  afforded  by  the  replacement  of  (105)  by  (106)  will  be  of 
interest  later;  see  especially  appendix  F,  (F-42)  et  seq.) 

Thus  the  leading  term  in  (95)  is  given  by  (105);  the  remaining 
terms  in  (95)  are  developed  below,  by  using  (74)  and  (69): 
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5,  sf  «*»)*  *  2.Ak  $  fe-fc*) 

sf  H*  »J  -  fix  ^0  ^ 

5^~Sy  rt,  nj  *  **M)  ^ 

3J1J  mT»?  *  I?!  (Vi)  ^N>  £(*»-*»),  (107) 


and  the  last  term  in  (95)  is  given  by  (86).  Thus  F  in  (95)  is  given  by  the 
sum  of  (105),  (107),  and  (86).  Employing  these  results  in  (93),  there  follows, 
by  use  of  (74),  (77),  (72),  and  (76), 


*  !)(|  +  2(k1  (r>+  T‘) +  ft)  r« 


(108) 


where  we  have  defined 


•*k 

di,  (H.  ^)h~ 


(109) 


and 


(110) 


Both  of  these  quantities  are  smaller  than  unity;  in  fact,  they  can 
be  much  smaller  than  unity.  For  signal  absent,  only  the  first  term  in  (108) 
remains  (see  (95)),  and  (108)  reduces  to  (88). 

When  we  combine  (75)  with  (108),  there  follows 


Tp,  denotes  a  p-^  fold  integral,  with  p  numerator  ^-functions  and  j,  denominator 
-functions.  If  ^  is  zero,  the  notation  is  simply  Tj. 
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wtoJM  -  (ij  z + 2  ({J  t3 + (if (X*-T>)  om 

The  first  term  in  (111)  is  identical  to  the  square  of  (89);  the  other  two  terms 
are  due  to  the  randomness  of  the  signal. 

Approximation  to  T» 

In  order  to  evaluate  J3  and  I#  of  (109)  and  (110),  we  need  to  make 
some  approximations.  These  approximations  will  be  developed  by  re-consideration 
of  T2,  and  carried  over  to  the  cases  for  Ta  and  I*, ,  where  exact  evaluation 
is  very  difficult.  First,  in  (76),  let  ■««  VlVOA  £>-  )i=  '>*■  tuw 

V* 

=  -=pr  du,  |r(u,^)r  =  JJdu,  du*  r{u,) r|i/,)|^(w,-u^|  f 

where  rectangular  weighting 

r|u)  o,  otW*i»«.  I  • 

Let  X  =  u.-u*  in  (112);  then 

Xx  “  JJdu,  dv  rK)  r(u,-x)|^f^)|  -  J  d*  ^ (*)  ty) 1  , 

where 

c)>x(x)  s  Jdw  r(u)r(u-x) 

is  the  autocorrelation  function  of  rectangular  weighting  (113).  This  auto¬ 
correlation  is  triangular  and  extends  over  the  range  (-T,  t).  On  the  other 
hand,  the  effective  width  of  |^lx)|a  is  approximately  W/' .  Thus  if 

TW,  l,  l/’&Or  has  decayed  to  zero  long  before  <fcjx)  has  changed 

significantly  from  its  origin  value.  This  allows  (114)  to  be  expressed 
approximately  as 

Xj.  *  <fjo)Jdx  M*)f  a  x  >f  T*4  »  (116) 


(114) 

(115) 


(112) 


(113) 
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1°  by 


In  order  to  evaluate  (116),  we  will  approximate  signal  correlation 
a  Gaussian  function.  We  choose 


fr )  a  (-  ?WS  T^)  fir  all  T. 


(117) 


The  choice  of  scale  factor  has  been  made  to  satisfy  condition  (80);  thus  W*  is 
the  statistical  bandwidth  of  the  received  signal.  In  appendix  E,  it  is  shown 
in  (E-28)- (E-33)  that  (117)  is  in  fact  a  good  approximation  to  the  true  signal 
correlation  when  the  rms  frequency  deviation  is  somewhat  larger  than  the 
equivalent  bandwidth  of  the  frequency-modulating  process.  Furthermore,  even 
when  this  condition  is  not  true,  Ws  in  (117)  could  be  chosen  so  that  (117)  fits 
fcs  well  as  possible  in  some  sense}  to  the  exact  signal  correlation,  given  by 
(E-ll)  and  ( E-9 ) .  Several  examples  for  modulating  and  signal  spectra  are  given 
at  the  end  of  appendix  E,  and  the  approximation  (117)  is  demonstrated. 

Substituting  (117)  in  (116),  we  find 

x-  s  fk  V  Tvt »  i ;  die) 

this  result  is  consistent  with  (79)  et  seq.  TWS  does  not  have  to  be  too  large 
for  (118)  to  be  a  reasonable  approximation. 

Equation  (118)  is  a  fairly  good  approximation  to  the  exact  curve 
for  c  =  0  in  figure  3,  where  the  signal  correlation  was  exponential;  see  (81). 

When  the  signal  correlation  f  is  Gaussian,  as  given  by  (117),  the  exact  value 
of  X,  is  available  from  (76)  according  to 


T  =  2^? 
ex' 


■[fW-i- 


I  -  exp(-*^ 


o<  =  y/Sr  TWS 


(119) 


Error  function  $  is  defined  in  appendix  E,  example  2.  Comparison  of  the 
approximate  and  exact  values  of  X*,  given  by  (118)  and  (119),  is  presented  in 
figure  4.  It  is  seen  that  (118)  is  a  rather  good  approximation,  especially  for 
large  TWS  .  Thus  the  approximation  technique  developed  in  (112)-(116)  will  be 
used  to  approximately  evaluate  T>  and  Tj*,  given  in  (109)  and  (110).  Actually, 
all  the  above  has  been  done  merely  to  set  the  stage  for  the  approximate  evaluation 
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of  X3  and  Tfc>  since  exact  evaluation  of  I,  via  (76)  is  rather  straightforward 
in  most  cases;  see,  e.g.,  (83)  and  (119).  However,  exact  evaluation  of  (109) 
and  (110)  can  be  extremely  difficult. 

Approximation  to  X» 

Use  of  (113)  allows  us  to  express  (109)  as 


Ij  =  r(u,)  ^Ux)r|«3)^fUarU,yilu,-u,)^?(u5-u^. 


where 


Letting  x  =  u»- u,,  ^  =  u»>  there  follows 

Ij  -  j|{  ')*  r|i)J  irlu,  hjy 

where 

={  cju  r(u-v)r(.u)>r(u+^). 

Now  if  T^»},  (12t>,  (122),  and  (113)  yield 

I3  *  fa(o,o)JJc)x  a(ty(-x-y)^(y)  =  ^rj \ ^  ^)jO(~^'^(y), 

since  the  integrand  of  (123)  — »  0  no  matter  how  x}y  — *  ±  <*5  ,  Next  we 
appeal  i»  approximation  (117)  in  order  to  evaluate  (123);  we  have  then 

Ti '  r ■ 

The  last  step  was  made  by  use  of  K-fold  multiple  integral  relation 
r  r  t  -»  /  K 


JdX  e*p[-XTMX]  - 


M 


,  dct  M  >  0, 


(120) 


(121) 


(122) 


(123) 


(124) 


(125) 
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where  column  matrix  X  *  [x,  ••  •  xK"]  ,  and  M  is  a  symmetric  K  x  K  matrix. 
Equation  (124)  constitutes  our  desired  approximation  to  X, . 


Approximation  to  X^ 

If  we  try  to  approximate  in  (110)  directly,  as  above,  we 
encounter  the  difficulty  that  the  appropriate  correlation  does  not  decay  to  zero 
for  large  arguments.  Specifically,  if  X, * and  in  (110),  the 

integrand  remains  at  value  1.  Furthermore  the  variance  expression  of  interest, 
(111),  involves  the  combination  Evaluation  of  this  latter  quantity 

is  undertaken  in  appendix  F;  the  end  result  )s 


(ntf 


(126) 


where  Xp  is  a  dimensionless  quantity  of  the  order  of  1 .  A  short  list  of 
values  of  X*  is  given  in  table  1,  for  Gaussian  frequency-modulating  process 
with  an  exponential  correlation  the  frequency-modulation  index 


D  = 


5l 


(127) 


is  the  ratio  of  the  RMS  frequency-deviation  to  the  equivalent  bandwidth  of  the 
frequency-modulating  process  vn(fc).  A  larger  table  and  plot  are  given  in  appendix  F 


0 

1  '5 

1 

2 

3 

4 

5 

6 

|  .356 

.418 

.642 

.881 

1  .124 

1  .368 

1.612 

Table  1.  Values  of  for  Exponential  Correlation 
of  Frequency  Modulating  Process  m(t) 


Approximation  to  Variance 

When  we  combine  (118),  (124),  and  (126)  in  (111),  we  obtain  the 

approximation 


(128A) 


where 


e/m.  _  aTz 


tw. 


M.W, 


_ received  signal  pouxar 

fzxeiv&i  t7*i5e  pouter iy  Jtjml  botnd 


(128B) 


is  the  input  signal-to-noise  power  ratio  (SNR)  in  the  signal  bandwidth*.  For 
very  large  TY/S  ,  keeping  E/l fixed  for  the  moment,  the  first  term  in  (128A) 
dominates;  that  is,  for  low  input  SNR,  the  first  term  is  the  same  as  the  square 
of  (89)  with  approximation  (118).  In  terms  of  deflection  4  defined  in  (90),  we 
can  employ  (92)  and  (118)  to  write  the  requirement  on  E//^  as 


■jE-  =  *  4  /rx  fc-  kg«.  TW,,  029) 

in  which  case  (128A)  and  (128B)  become 

Vorfttf  W)]  =  ^  [' +  ^  +  T>  J  (130A) 

and 

_  e/m,  _  j, 

‘  T*k  "  /Rf  '  (,30B) 


Another  interpretation  of  R,  is  furnished  by  expressing  ??*®E,/A/, ,  where  E,*  tvif 
is  the  signal  energy  per  independent  component.  That  is,  in  observation  time  T, 
there  are  tW4  effectively  independent  components  of  the  received  signal  (see 
the  footnote  to  equation  (80)).  Then  small  "Rs  is  synonymous  with  a  small  ratio 
EJbi,  of  the  component  signal  energy  to  noise  density  ratio.  This  interpretation 
will  be  of  use  later. 
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The  first  term  in  (13QA)  is  the  only  one  that  would  remain  for  hypothesis  H,, 
signal  absent;  see  (89).  Thus  the  increase  caused  by  signal  randomness  is 
indicated  by  the  last  two  terms  in  (128A)  or  (130A),  in  terms  of  input  measure  % 
or  output  measure  <i»,  respectively.  If  yFfw*  is  somewhat  larger  than  the 
performance  quality  measure  d*  (see  (90)- (91)),  then  (130A)  is  only  slightly 
larger  than  the  variance  with  signal  absent. 

Although  the  variability  of  £")£*  is  somewhat  greater  with  signal- 
present  than  with  signal -absent,  this  may  not  matter  much  if  we  are  asking  for 
detection  probabilities  in  the  neighborhood  of  0.5.  For  this  case,  the  threshold 
will  lie  approximately  in  the  center  of  the  probability  density  function  of  the 
decision  variable  when  signal  is  present,  and  the  increased  variability  is  not 
relevant.  And  for  false  alarm  probabilities  in  the  )0  range,  with  TYt'*), 
the  deflection  criterion  of  (90)  is  the  appropriate  quantity  to  focus  on.  Thus, 
if  we  confine  attention  to  the  i:  ^  term  of  the  In  LR  development  in  (49), 
requirement  (129)  is  the  one  of  major  significance;  whether  higher-order  terms 
are  important  has  not  yet  been  ascertained. 


DERIVATION  OF  *r 

We  are  now  prepared  and  interested  in  determining  some  of  the 
higher-order  terms  in  the  series  expansion  of  In  LR  in  (49).  We  recall  that 
r  is  given  by  (32)  or  (42),  and  that  and  yf?  are  given  by  (39). 

Substituting  (96)  into  (32),  we  have 

r  =  1?«{Ce*p (-»>)}  s  cVf),  (131) 

where  we  have  defined  the  complex  quantity 

(2  =  -n-Jdt  X  (0  w{u)j  (132) 

Now,  according  to  (39),  we  must  average  r3  over  the  signal  statistics,  i.e., 
over  ^  and  m(t)  Performing  the  average  over  ^  first,  we  obtain  from  (131) 

(rVK(^^V»y>ro,  0  33, 

since  <|>  is  uniformly  distributed.  Therefore 

V3r>  =  0,  (134) 

meaning  that  the  term  in  LR  of  (49)  is  absent.  In  fact,  all  cumulants, 

f  for  k  odd,  are  zero,  by  an  argument  similar  to  (1 31 )- (1 33) ,  for  the  signal 
model  (96). 
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DERIVATION  OF  ^ 

We  now  have  to  average  r+  over  the  signal  statistics.  From 


(131),  we  find 


’|cr 


(135) 


Hence,  from  (39), 


A<icr>„  -^xr 


(136) 


here  «  must  still  average  Icf  over  the  random  signal  frequency-modulating 
process  >*K)-  We  have,  from  (132), 

Jum(u)- i  • 


(137) 


But  we  have  already  evaluated  this  statistical  average,  in  (97)-(105).  Therefore 


we  find 


<Kf >.  -  ©  JflK***  • 


(138) 


where  we  have  also  used  (77). 

Combining  (138),  (45),  (74),  and  (75)  in  (136),  we  obtain 

X--f  =  -j-  /At fa  r  &  dt.  Jt,  Jt  7S.lt)  A)  K  U  2*4,)  • 


(139) 
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This  is  the  next  non-zero  term  in  expansion  (49).  Although  it  Indicates 
explicitly  what  has  to  be  done  on  the  received  waveform  no  simple  way 

of  realizing  this  processing  has  been  discovered;  in  general,  it  appears  that  the 
difficult  processing  indicated  by  the  fourth-order  integral  must  be  evaluated. 
Whether  there  is  a  special  property  of  the  combination  of  correlations  in  (139), 
that  enables  a  simpler  or  approximate  realization,  is  unknown.  The  magnitude 
of  the  difficulty  of  calculating  (139)  directly  is  pointed  out  by  noting  that 
since  samples  must  be  taken  in  time  at  least  as  often  as  seconds,  and  xW 

is  of  duration  T,  a  sum  of  terms  of  order  must  be  effected.  For  large  TV,, 
this  is  not  going  to  be  possible,  in  which  case  approximations  must  be  found. 

For  now,  however,  we  are  interested  in  the  average  size  of  the  term  (139)  in 
the  In  LR  series. 


A  partial  check  on  (139)  is  afforded  by  using  the  example  of  xH)  = 
which  leads  to  *4 /24  =  -(V'V/H.  At  the  same  time,  (132)  then  yields 

ICl*  =  AVMo  ,  while  (131 )  yi el ds  <r*>  =  f aVa£  and  <r4> -- and 

(39)  yields  the  same  value  as  above  for  7^/2+.  Alternatively,  use  of  the 
replacement-approximation  of  (105)- (106)  results  in  zero  for  (139);  this  is 
consistent  with  a  Gaussian  assumption  for  sH).  More  precisely,  use  of  (F-47) 
in  (139)  yields 


where  remainder  R3  tends  to  zero  no  matter  how  increase  to 

infinity.  ' 
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MEAN  VALUES  OF 


Signal  Absent 

Let  the  combination  of  correlation  functions  in  (139)  be  denoted 
by  kernel  K  ±,,-{4)  (which  equals  -H)  V  ^*)  )  •  Then 

for  noise-alone,  using  (85)  and  (86),  along  with  (76),  we  have  mean  value 

,4 


I 

14- 


Jif{  *<<**>  4  ^  • 

[R-  t»l  $  Ry  +,)  +  Ur  L)  l 


-  .i£ 
~  t 


^  II  ^  ^  Mj)  I  =  i  (mJ  Ij 


(140) 


For  comparison,  (71)  yields  the  mean  value  of  the  previous  term 
in  the  \n  LR  expansion  (49): 


(141) 


Since  T,.  r  (tW)  for  via  (118),  (140)  becomes  small  relative  to 

(141),  as  TV,  increases,  as  was  anticipated,  provided  that  E/^e  is  held 
constant;  however,  large  values  of  TW,  would  be  required,  due  to  the  additional 
power  of  E/M,  involved  in  (140).  Furthermore,  the  situation  of  fixed  £/N„  is 
not  the  one  of  interest;  rather,  if  we  attempt  to  keep  the  deflection  criterion 
dj  of  (90)  constant,  we  have  from  (92), 

xT(n)  *  -  i  ^  •  <142> 


which  is  constant  (and  large),  while 


*  J./m'  V  TW,»/. 


(143) 


In  order  that  (142)  be  much  smaller  (in  magnitude)  than  (143),  we  would 
require  HA/twT  =i*s  «  I,  which  would  indicate  large  values  for  TW,, 
since  values  of  Ji/io  are  typical  requirements  for  good  detection 
performance;  see  the  discussion  surrounding  (90)- (91 ) .  However,  these  mean 
values  are  not  the  final  or  most  meaningful  measures  of  performance;  additional 
important  parameters  are  evaluated  below. 

Signal  Present  ^ 

The  derivation  of  the  mean  value  of  2<f  in  (139),  for  signal 
present,  is  rather  involved  and  is  presented  in  appendix  G;  the  result  is 

feT  k ' 2  T‘^ 


(144) 


where 


Mr-Qvk' 


(145) 


and  T,,r,,  and  are  given  by  (76),  (109),  and  (110).  Evaluation  of  1^-23^ 
in  appendix  G,  and  use  of  (118)  and  (124),  yields 


*555 gr  ^ 

+  '  i1*  Rsl]  ^  ’“T  TW” 


(146) 


(147) 


-rw 

where  is  a  dimensionless  quantity  of  the  order  of  unity.  In  terms  of 

deflection  d2,  we  use  (129)  to  express  (146)  as 
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Comparison  of  (142)  and  (148)  reveals  that  the  difference  of  mean  fourth-order 
outputs  is 


.  <’49> 

which  can  be  large  for  /io’-  However,  when  we  recall  that  the  difference 

of  mean  second-order  outputs  is 

J- M  -  - (0 *,-£>  nso) 

where  we  used  (78)  and  (92),  then  we  see  that  the  latter  difference  is  much 
larger  if  1?s  =  is  small.  But  this  is  still  not  the  final  statement, 

for  we  still  have  to  address  the  deflection  criterion  of  random  variable  3^, 
which  will  involve  the  variance  of  ^(w)  ,  in  addition  to  the  difference 
of  means  considered  above. 


i 

i 
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VARIANCE  OF  ^(n) 

Since  ^  is  given  by  the  four-fold  integral  in  (139),  evaluation 
of  its  variance  will  involve  an  eight-fold  integral.  To  simplify  this  task,  we 
will  consider  the  variance  only  for  noise-alone.  The  derivation  Is  given  in 
appendix  H,  with  the  result  that 

-  i(£)  (r«-  +< 3;+2x;),  on) 


where  the  various  quantities  have  been  defined  in  (76),  (110),  (145),  and  also 


4,  dt.4  A,  A-k)  f  i) 


(152) 


The  expression  (151)  is  simplified  in  appendix  H,  to 


(153) 


_tfl 

where  i-p  is  a  dimensionless  quantity  of  the  order  of  unity,  and  is  available 
from  ijj*  and  according  to  (H-8).  In  terms  of  deflection  d*,  (129)  yields 


4- 


(154) 


For  comparison,  the  variance  of  the  first  term  in  the  )n  LR  series 
is  given  by  (89)  and  (118)  as 

Varfi'VfW}  =(;f-)  ^  ■  0  55) 

The  ratio  of  (154)  to  (155)  is  of  the  order  of  ( d,  /  -  Hi  ,  which  will  be 
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rather  small  If  /Tty*  is  somewhat  larger  than  <1*. 

We  can  now  form  a  deflection  statistic  for  the  second  term  alone, 
r)^/24  ,  In  the  In  LR  series.  Using  (140),  (146),  and  (153),  it  is 


i&W 


.  056) 


In  terms  of  deflection  we  use  (129)  to  find 


(157) 


Thus  the  second  term  in  the  In  LR  series  has  about  the  same  deflection  criterion 
as  the  first  term,  when  fn£  is  large  relative  to  deflection  .  At  first, 
this  result  would  seem  to  suggest  that  the  In  LR  series  ought  to  be  carried  to 
at  least  two  terms;  however,  a  better  measure  for  this  conclusion  is  based  upon 
the  deflection  criterion  of  the  combined  statistic 


(158) 


rather  than  upon  each  term  alone.  We  now  undertake  the  evaluation  of  the 
deflection  of  statistic  (158). 
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DEFLECTION  OF  J0* 


The  difference  of  means  of  for  signal  present  versus  noise-alone, 
can  be  obtained  from  (78),  (118),  (140),  and  (146);  we  find,  with  the  help 


of  (1308), 


(159) 


The  variance  of  for  noise-alone  is  given  by 

+  2  .  (160) 


With  the  exception  of  the  third  term,  these  quantities  are  available  from  (89), 
(153),  (71),  and  (140).  The  third  term  is  evaluated  in  appendix  I,  with  the 
result 

A(n)jUn)  -  =  ~(j|)  X3‘  (161) 


Combining  this  result  with  those  cited  above,  and  employing  (118),  (124),  and 
(130B),  there  follows 


Vor^M]  =  (f)  f>4 


■h 


-1(£ f 

+  VNo  / 


(162) 
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The  deflection  of  Jj*  is 


SD  [ij*')') 


e_  _j _ 

^  ^  +-oU0' 


where  we  employed  (159)  and  (162),  and  recall  (130B), 


E _ [_  _  cL 

Me  T\WS  i/t>£ 


Expanding  (163)  in  terms  of  powers  of  T?s,  we  find 


_w[i  +  oW]  .4*0(4.)} 


(163) 


(164) 


(165) 


Notice  that  the  perturbation  from  is  an  order  of  magnitude  smaller, 
depending  on  (d, /\Jt*Js  ')  now.  Thus  the  addition  of  ^  to  A  does  not 
significantly  affect  the  deflection  criterion;  this  strongly  suggests  that  the 
first  non-zero  data-dependent  term  in  the  In  LR  series,  Jlx  =  "t  hO? )  is 
sufficient  to  base  optimum  decisions  on,  without  regard  for  the  higher-order 
terms,  when  {t\  >  3 and  di^v/To*.  This  criterion  is  significantly  weaker 
than  the  earlier  ones  which  were  based  on  preliminary  statistics,  and  not  a 
final  performance  measure. 

Correlation  Coefficient 

One  additional  statistic  relating  Jl  and  is  of  interest, 

* 

namely  their  correlation  coefficient;  for  noise-alone,  it  is  given  by 


i - - 

Calculation  of  the  correlation  coefficient  for  signal-plus-noise  would  require 
knowledge  of  the  mean  of  J1,(S*M)  -fA(5f  n)  and  the  variance  of  ^(Sf  n),  in 
addition  to  the  other  quantities  that  have  been  evaluated  already.  This  would 
mean  evaluation  of  sixth-order  and  eighth-order  signal  and  noise  correlations  and 
subsequent  simplification  and  approximation;  this  task  has  not  been  undertaken. 
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,  ,  _  frl)  -  A  frj) 

~  sd&Iu)]  s»  UM] 

Substitution  of  (161),  (89),  (153),  (118),  and  (124)  yields 


(166) 


■for  )ci»^c.  TWs  . 


Evaluation  of  Tp  is  undertaken  in  appendix  J;  we  find  via  (J-13) 


(167) 


C  Cm)  -  -  W  b  =  3-  dee) 

This  rather  large  correlation  coefficient  indicates  that  the  statistics  J2(n) 
and  intend  to  vary  together;  that  is,  when  one  gets  large,  so  does  the 
other  (in  magnitude),  and  vice  versa.  This  result,  combined  with  (165)  above, 
is  further  confirmation  that  the  addition  of  £+  to  does  not  significantly 
alter  the  decisions  yielded  by  £  alone.  This  conclusion  is  drawn  for  large 
TWS  only.  As  an  example,  if  we  desire  \f\o  (see  (90 )- (91 ) ) ,  then  (165) 
suggests  T>JS>IOO  as  the  condition  for  neglect  of  .  Consideration  of 
(130B)  reveals  that  we  are  then  talking  about  input  SNRs  less  than  //✓To’  = -5 'JR 

To  summarize,  if  we  want  decent  performance  in  terms  of  detection 
and  false  alarm  probabilities  (d1£v/7o'1),  and  if  /rWs'  is  large  relative 
to  d,  (yTVds’  k  ,  then  the  input  SNR  is  small  (1?s  <  '/j),  the  time-bandwidth 
product  is  large  (tW5 >  / oo),  and  the  4-  term  1n  the  In  LR  series  can  be 
dropped  (use  fLz  alone).  In  fact,  if  ^rw5'  >  Max  ( lo,  3  ,  but  d,  is  small, 

then  'R*<  •/>  >  an^  the  term  can  still  be  dropped  without  degradation  in 
performance;  however  the  performance  level,  in  terms  of  false-alarm  and  detection 
probabilities,  will  be  poor. 

If  we  attempt  to  state  the  condition  for  neglect  of  in  terms  of 
input  SNR  1?s,  we  can  say  that  if  1?s  <  »/3,  then  from  (164),  fvNg  >3  da .  However, 
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this  does  not  necessarily  mean  that  TV^  is  large  unless  3.  So  small 

input  SNR  is  not  a  sufficient  condition  by  itself;  it  must  be  augmented  with 
the  requirement  of  good  performance  by  the  processor.  Then  we  can  state  that 
can  be  dropped  in  the  LR  series. 


ON  THE  APPROACH  OF  r  TO  GAUSSIAN 


It  was  noted  in  (47)  et  seq.  that  the  random  variable  r  should 
approach  Gaussian  as  TWS  increases;  we  can  now  make  some  quantitative  statements 
about  the  rate  of  approach.  If  r  were  equal  to  a  sum  of  N  independent  identically 
distributed  random  variables  [%.},  ,  then  we  would  find  that  the  measure  of  non- 
Gaussianness  afforded  by 


2?  l 

if  w  X*' 


(169) 


indicates  a  '/n  decay.  For  our  problem,  however,  r  is  given  by  (47)  as  a  sum 
of  approximately  TWS  unequally-distributed  independent  components;  the  unequal 
distribution  follows  from  the  weighting  by  the  received  waveform  x(D,  which 
varies  with  time  t\  Furthermore,  this  results  in  themselves  being  random 

variables.  Nevertheless,  we  can  still  derive  a  meaningful  measure,  similar  to 
(169),  for  our  problem. 

We  observe  first  that  the  random  variables 
are  clustered  if  I  ;  that  is,  their  standard  deviations  are  much  smaller 

than  their  means,  if  TWS»|.  From  (71),  (89),  and  (79), 


l/m 


(170) 


From  (75),  (79),  and  (128), 


where  is  defined  in  (164);  the  function  of  Rs  is  approximately  unity.  And 
from  (140),  (79),  and  (153), 
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072) 


the  numerator  of  (172)  is  approximately  3.44  for  D  =  3;  see  (J-13).  For 
large  TW5 ,  these  last  three  quantities  are  small,  and  our  conclusion  about 
clustered  behavior  is  drawn. 

This  leads  us  to  replace  (169)  by  the  measures 


jffjjo 

M 


(173) 


(174) 


which  are  obtainable  respectively  from  (140),  (79),  (71),  and  (146),  (75), 

(79),  (164).  Since  the  function  of  l?s  in  (174)  is  approximately  unity,  both 
measures  (173)  and  (174)  indicate  a  rate  of  decay  proportional  to  (TH)  J 
this  inverse  dependence  on  the  effective  number  of  independent  components  in  r 
is  similar  to  (169).  The  scale  factor  in  (173)  and  (174)  is  somewhat  larger 
and  is  probably  due  to  the  unequal  weighting  in  (47),  as  mentioned  above. 

The  cumulants  ^  and  ^  are  zero,  whether  signal  is  present 
or  not,  as  are  all  odd  cumulants-,  thus  is  the  first  non-Gaussian  contributor 

to  the  In  LR  series: 


UK- 


-7f +••• 


(175) 
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where  we  employed  (49),  (41)  et  seq.,  and  (134).  All  of  the  above  measures 
lead  to  the  conclusion  that  the  most  important  term  in  (175)  is  the  first 
data-dependent  term  ,  at  least  for  large  TWS. 


POWER  AND  ENERGY  RELATIONS 


1 


It  has  been  noted  above  that,  for  large  TW, ,  the  dominant  term 
of  the  U  LR  series  is  the  term,  and  that  in  order  for  this  term  to 

realize  deflection  d2)  we  must  have  the  ratio  of  received  signal  energy 
to  noise-density 


(176) 


where  ^ ~  /fi?  for  good  performance;  see  (129).  This  relation  requires  the 
received  signal  energy  E  to  increase  according  to  the  square  root  of  the  signal 
bandwidth,  Ws  ,  for  fixed  deflection  noise  density  A4 ,  and  observation 
time  T- 

The  ratio  of  the  received  signal  power  to  the  noise  power  in  the 
signal  band  is  then 


,  _  K/z  _  e/t  _  dg 
s  H  Ws  /-rw7 


(177) 


which  allows  to  decrease  as  Ws  for  fixed  d*  and  T.  The  reason  for 
this  decrease  is  that  the  received  noise  power  in  the  signal  band  is  increasing 
according  to  Ws,  as  Ws  increases.  A  better  measure  may  be  the  ratio  of 
received  signal  power  to  the  noise  power  in  a  1  Hz  band: 


Like  E/N,,  this  must  also  increase  as  W ,  for  fixed  and  T.  Some 
examples  of  these  relations  are  presented  in  table  2. 


^a 

w5W 

TH 

TWS 

■RsW 

VTo-1 

31 

*3 

300 

1000 

20 

-10 

-4.8 

fuT 

1 

3 

300 

100 

15 

-  5 

-9.8 

lfT<r 

1 

30 

300 

* 

10 

10 

0 

-14.8 

Table  2. 

Power  and 

Energy  Relations 

The  input  signal-to-noise  ratio,  (177),  can  range  below  0  dB  if 
large  values  of  TW5  are  attainable;  but  rather  large  values  for  E/At  are  required 
in  this  case.  For  comparison  and  partial  verification  of  the  last  example  in  table 
2,  Ref.  8,  page  14  gives,  for  M  ~  1  alternative^  =  10-fold  diversity,  ?p  =  10  *  £  =. S', 
the  value  Er/A4  =  11.3  c)B;  thus,  high  orders  of  "diversity"  (TW,*>  l)  require 
large  values  for  the  received  signal  energy  to  noise  density  ratio,  in  order  to 
attain  decent  performance. 


This  is  not  a  very  large  value  of  TWS  >  thus  the  conclusion  about  = 
being  a  sufficient  statistic,  and  the  adequacy  of  as  a  performance  measure, 
are  suspect  for  this  last  example. 
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THIRD  CUMULANT  OF 

The  first  two  cumulants  of  have  been  derived  in  (71) 

and  (89);  they  are,  using  (90), 


-  c. 


=  (0  Tl  =  ^  ’ 


(179) 


For  large  TW5(  these  relations  were  the  basis  in  (170)  for  showing  that 
X  is  clustered  around  its  mean  value.  Now  we  wish  to  show  that 

i:  In)  is  tending  to  a  Gaussian  random  variable  as  TWS  increases.  (This 

is  true  even  though  ^  X^In)  is  a  nonlinear  double  integral  of  the  noise 
waveform.)  We  do  this  by  considering  the  normalized  third  cumulant 
which  is  a  standard  measure  of  non-Gaussianness  of  a  random  variable  (see,  e.g. 
Ref.  4).  We  have 


C3  -  [tVTW- 


(180) 


We  expand  (180)  out  and  use  (71)  and  (88)  to  evaluate  the  first-  and  second-power 
moments.  The  third  moment  can  be  evaluated  by  using  the  properties  of  complex 
white  Gaussian  noise  in  (45),  with  the  result  that 

^3  =  2  >  (181) 

where  X,  is  defined  in  (109).  Therefore  (179),  in  conjunction  with  (118)  and 
(124),  yields 


=  2  -5>  w  TWS  (182) 

l/‘  i/twT  17 

The  square-root  decay  with  the  number  of  independent  contributors  is  again 
standard  for  this  particular  normalized  cumulant;  thus  i  is  tending 

towards  Gaussian.  This  result  (182)  could  also  be  useful  in  setting  thresholds 
for  specified  false  alarm  probabilities  of  the  processor  employing  only 
in  its  decision-making. 
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AND  ITS  MEAN 


SIXTH  CUMULANT  ^ 


The  conclusions  drawn  above  have  been  based  on  consideration 
of  and  and  their  statistics.  We  now  consider  the  sixth  cumulant 

y*1  and  its  mean  ^(n)  with  noise-only;  any  other  statistics  require  a 
prohibitive  amount  of  manipulations.  The  derivation  is  presented  in  appendix 
K,  with  the  result 


X.  Vr) 


23  04 


[*&> -A) -Kfc,  A)  KM)  +  UK.  MMM.M)]  > 


(183) 


where 


K,  0  S  <f*lt ,  (184) 


is  the  complex  phasor  process  defined  in  (F-43),  is  available  in 
(F-44) ,  and  Ki=/°-  The  complex  processing  required  by  the  kernel  of  (183)  has 
not  been  evaluated,  because  it  is  too  difficult  to  realize  physically  and  sheds 
no  light  on  the  processor. 

The  mean  of  (183)  for  noise-alone  is  derived  in  appendix  K: 

■  i®*. '  f-w™. 


(185) 


For  comparison,  the  means  of  the  first  three  non-zero  data-dependent  terms  in 
the  in  LR  series  are  summarized: 


B  4  p* 
4  3&  * 


(186) 


Thus  the  mean  values  of  the  higher-order  terms  are  dropping  off  by  additional 
powers  of  the  input  SNR  which  is  small.  This  leads  to  the  conjecture  that 
the  higher-order  terms  in  the  In  LR  series  are  progressively  less  important  when 
w  >  3  ,  i  .e. .  small  7?,. 


SUMMARY 


The  logarithm  of  the  likelihood  ratio  has  been  developed  in  an 
infinite  Vol terra  series,  and  the  leading  terms,  from  zeroth-order  through 
sixth-order,  have  been  investigated  in  terms  of  their  size  and  statistical 
significance.  It  has  been  found  for  the  fully-random  signal,  which  is 
characterized  by  no  deterministic  component  and  no  carrier,  that  if  false 
alarm  probabilities  in  the  order  of  )0*3  and  detection  probabilities  in  the 
order  of  .5  are  desired,  if  the  product,  TWf,  of  observation  time  and  received 
signal  bandwidth  is  large,  and  if  there  is  at  least  a  moderate  frequency  modulation 
index  D,  then  the  ratio  of  total  received  signal  energy  to  noise  power  density 
level  must  be  rather  large,  of  the  order  of  3/fW?  ,and  the  dominant  term  in  the 
log  likelihood  series  is  given  by  the  filter-energy-detector  term,  The 

requirement  on  the  frequency  modulation  index  D  being  at  least  of  the  order  of 
2-3  was  introduced  in  the  analysis  in  order  to  facilitate  approximate  evaluation 
of  some  of  the  multiple  integrals  by  means  of  Gaussian  functions.  Thus  this  is 
a  sufficient  condition  employed  for  tractability;  how  small  the  modulation  index 
can  become,  without  violating  the  conclusions  above,  is  unknown. 

Alternatively,  if  VfwT  is  larger  than  the  maximum  of  10  and  3di, 
where  4*  is  the  deflection  criterion  of  the  filter-energy-detector  term  ^ , 
then  the  dominant  term  in  the  log  likelihood  series  is  }£’ ;  this  holds  regardless 
o^  the  size  of  42,  i.e.,  whether  good  or  poor  performance  is  to  be  obtained  from 
the  filter-energy-detector.  In  this  case,  the  input  signal-to-noise  power  ratio 
in  the  signal  band  is  small.  It  is  not  sufficient  to  say  that  small  input  signal- 
to-noise  power  ratio  alone  yields  the  filter-energy-detector  as  near-optimum. 

Rather  this  condition  must  be  augmented  with  a  requirement  for  large  TW,,  i.e., 
/T>£  >(l©,3<U),  which  means  that  larger  time-bandwidth  products  are  required 
for  better  performance. 

The  input  signal-to-noise  ratio,  can  be  put  in  the  form 
1?5=  E,/hl»,  where  E|  is  the  received  signal  energy  per  independent  signal 
component  and  is  the  noise  power  density  level;  see  the  footnote  to  (128B). 

When  this  ratio  of  component  signal  energy  to  noise  density  is  small,  and  the 
number  of  components,  7WS,  is  large,  the  optimum  processor  is  well  approximated 
by  the  filter-energy-detector.  This  same  conclusion  has  been  reached  by  the  author 
in  some  as  yet  unpublished  work  (Ref.  9)  on  the  exact  performance  of  a  related 
processor,  the  combiner;  see  (51)  et  seq.  Specifically,  for  a  signal  with  K 

orthogonal  signal  components,  deterministic  except  for  phase,  the  optimum  processor 
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m 


is  a  sum  of  In  of  the  sampled  envelopes  of  the  outputs  of  K  filters,  each 
one  matched  to  a  different  signal  component.  When  the  number  of  Independent 
components,  K,  tends  to  infinity,  the  optimum  processor  does  not  approach  a 
sum  of  envelopes-squared  (the  filter-energy-detector)  unless  the  signal  energy 
to  noise  density  ratio  of  each  component  tends  to  zero  as  K  increases.  Also, 
in  this  case,  the  performance  of  the  optimum  processor  and  the  sum  of  envelopes- 
squared  processor  are  close  to  each  other  in  the  region  of  reasonably-good 
false  alarm  and  detection  probabilities.  This  exact  back-up  analysis  of  a 
related  processor  tends  to  confirm  the  results  achieved  in  the  current  study. 

We  have  assumed  here  that  the  signal  amplitude.  A,  is  a  constant 
over  the  observation  interval.  This  led  to  a  more  tractable  analysis  and  did 
not  require  any  additional  assumptions  about  amplitude  statistics  and  their 
dependence  on  frequency-  or  phase-modulation  statistics.  It  is  this  author's 
conjecture  that  in  the  case  of  random  amplitude-modulation  in  addition,  the 
filter-energy-detector  will  again  turn  out  to  be  near-optimum  under  the  conditions 
of  small  component  signal  energy  to  noise  density  and  a  large  number  of  independent 
components.  This  is  based,  in  part,  on  the  observation  that  the  received  signal 
would  be  even  more  random  than  in  this  study,  and  that  the  filter-energy-detector 
is  a  robust  processor  for  the  more-noise-like  signals.  Of  course,  the  pre-filter 
in  the  energy  detector  (see  figure  2)  would  have  to  be  broadened  to  cover  the  total 
bandwidth  of  the  received  signal  (with  both  AM  and  FM),  thereby  lowering  the  input 
signal-to-noise  ratio  in  the  received  signal  band.  This  would  cause  a  loss  in 
performance,  but  is  unavoidable  as  the  randomness  of  the  received  signal  is 
increased.  Alternatively,  to  maintain  a  desired  level  of  performance,  the 
observation  time  would  have  to  be  increased;  the  result  is  a  larger  number  of 
independent  signal  components  and  a  still-better  approximation  to  the  optimum 
processor. 

Due  to  mathematical  difficulty,  it  has  not  been  possible  to  develop 
anything  but  an  infinite  series  for  the  log  likelihood  ratio.  We  then  had  to 
analyze  the  low-order  terms  (through  order  seven  which  is  zero)  to  determine 
which  were  significant,  and  under  what  conditions  the  leading  data-dependent 
term  was  dominant.  Every  statistic  (we  could  reasonably  evaluate)  pointed  to  the 
filter-energy-detector  as  being  the  dominant  term,  under  the  conditions  cited 
above.  Evaluation  of  higher-order  statistics  of  XiT*  such  as  their  higher- 
order  cumulants  (see  (182),  for  example)  or  higher-order  terms  (i.e.,  larger  k, 
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as  In  (183)  and  (186),  for  example)  Is  possible,  but  Is  extremely  tedious.  It 
also  requires  approximations  to  still  higher-order  multiple  integrals,  which 
have  proven  very  time-consuming  to  evaluate  exactly.  Furthermore,  the  difficulty 
of  realizing  Volterra  kernals  such  as  (139)  and  (183)  appears  to  preclude  this 
series  as  a  practical  solution  to  the  problem  of  optimum  detection  of  frequency- 
modulated  tones.  Rather  it  appears  to  this  author  that  for  moderate  T\*4  products, 
approximations  to  the  likelihood  ratio  directly  should  be  attempted.  Analysis  of 
their  performance  will  probably  require  extensive  simulation  effort. 
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APPENDIX  A.  OPTIMUM  PROCESSOR  FOR  SMALL  TW, 


If  I,  we  can  give  an  explicit  expression  for  the  LR. 

From  (32), 

r  S  Q*]  kr  f  «  T  <  ^  *  (A— 1 ) 

where  2.  and  a  are  the  values  of  *H)  and  A 19  at  the  endpoint  V  Now  let 
us  consider  s  =  ls|  where  ^  is  uniformly  distributed  over  2-rr. 

From  (13),  since  E  •=  ^Ulx"n  then  ls|  can  not  be  random,  since  E  is 
fixed.  Now  we  have  deterministic  component  <s>  -o,  and  therefore,  from  (33), 


a  *  ^  *  Ul 


(A-2) 


Also  <v>=0  from  (31).  Then  (35)  yields  (since  is  uniformly  distributed) 

U  =  «*}])> 

=  ^(-f)T.(tls|ui)  <A-3> 

Hence  comparison  of  |x|  with  a  threshold  is  optimum. 

On  the  other  hand,  the  quantity 

£•)£-  (A-4) 


This  again  yields  the  rule:  compare  Ixl  with  a  threshold. 

This  rule  for  decision-making  is  physically  reasonable,  since  a 
single  sample,  of  a  complex  envelope  with  uniform  phase,  contains  only  one  item 
of  information  on  which  to  base  decisions,  namely  the  magnitude  of  the  complex 
envelope. 


APPENDIX  B.  SOME  SIGNAL  MODELS  AND  THEIR  INTERPRETATION 
Consider  the  phase-modulated  signal  process  (also  given  in  (14)) 

sftS  -  A  «s[arf.t  +  6ft)  +  <■<  ,  (B-1) 

where  phase-modulating  process  ett)  is  a  zero-mean  real  stationary  Gaussian 
random  process,  and  is  uniformly  distributed  over  2ir.  Then  complex  envelope 

yft)  =  A  (b-2) 

Obviously  hence  there  is  no  deterministic  signal  component.  The 

ac  signal  component  is,  from  (33), 

A  V  <B-3> 

Then  the  averages  over  j), 

=  o  -  -  o,  (b-4) 

as  expected.  Also 

<^a  ft,)  fl*  fe))>  -  A3  (  e*j>  [i  e  ft,)  - «  0  ft.)]  y& 

=  A  e*j>[-Ke(o)-t-  {B_5) 

where  we  have  used  the  Gaussian  character  of  the  difference  variable  6^0-0,10, 
and  defined  the  autocorrelation  of  stationary  process  6  ft)  as 

<9ft>)&ft,)^  =  Rek'-t>)-  (B-6) 

v 

So  we  have,  from  (B-5) 

K3 lo>  =  a’,  K,W'  A2 vxr [-)?,(!>)],  (B-7) 
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since 


K,M  =  o. 

«,(»)  *  <9'«>, 

is  large,  i .e. , 


The  latter  quantity  in  (B-7)  is  not  zero  unless 
is  infinite.  However  if  the  mean-square  phase-modulation 

<?«£>  l, 


then  (B-7)  yields 


2*0. 


(B-8) 


(B-9) 


But,  in  general,  a  small  mean-square  phase-modulation  will  give  rise  to  a  non¬ 
zero  value  for  This  holds  even  though  there  is  no  deterministic  signal 
component,  i  .e. ,  <” slt)^  =  0. 

For  the  case  of  frequency-modulation  rather  than  phase-modulation. 


we  have 


(B- 10) 


where  mH)  is  the  instantaneous  frequency-modulating  process,  assumed  stationary. 
This  form  allows  6$  to  be  a  non-stationary  process.  The  properties  of  Gaussian 
frequency-modulation  are  taken  up  in  appendix  E;  it  is  shown  that  the  exponential 
quantity  in  (B-5)  is  replaced  according  to 

=»  ,  (B-n) 


where  6rw({)  is  the  power  density  spectrum  of  stationary  frequency-modulating 
process  •»(+).  Hence  (B-5)  yields 

The  behavior  of  (B-ll),  as  T-+<w,  depends  on  the  behavior  of  W 
near  -f=0i  numerous  examples  are  given  in  appendix  E.  Suffice  it  to  say,  for 
now,  that  examples  exist  where  (B-ll)  tends  to  infinity,  and  other  examples  where 
(B-ll)  tends  to  a  finite  constant.  Thus  (B-12)  can  tend  to  zero  or  non-zero  values, 
respectively,  as  T-**.  Furthermore,  this  holds  even  though  the  mean-square 
frequency  modulation. 


(B-l 3) 


is  finite.  Thus  detection  of  frequency-modulated  processes  will  be  different 
for  the  two  cases  of  (B-7T )  being  finite  or  infinite  in  the  limit  as  T-*  °* ; 
see  figure  2  in  the  main  text. 

When  the  signal  has  no  deterministic  component,  i.e., 

*  0, 

and  in  addition 

M*)--  o, 

we  call  the  signal  fully-random;  this  is  the  case  of  major  interest  here. 
However,  below,  we  delineate  the  components  of  the  signal  in  the  general  case, 
so  that  we  can  properly  interpret  our  case  of  interest. 

Interpretation  of  Signal  Components 

The  narrowband  signal  is  represented  in  terms  of  its  complex 
envelope  according  to 

s-ft)  *  Re  {  s  ft)  exj>(i  2rr£-0]  ,  <  4  <  V 

The  ac  component  of  the  complex  envelope  is 

2  (£)  -  OWX 

Hence  the  signal  can  be  represented  as 

s|t)  -  (i2»t,-t)} 

=  dlt)  +  alt), 

where  d(t)  is  a  deterministic  narrowband  waveform  with  known  phase-modulation, 
centered  at  frequency  -f#.  This  component  could  include  pure-tone  components  of 
known  phase,  including  as  a  special  case,  a  carrier  at  -P, ,  and/or  sidebands  in 
the  neighborhood  of  P0. 


(B-14) 

(B-15) 

(B-16) 

(B-17) 

(B-18) 
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The  correlation  of  signal  s|0  is,  using  (B-18), 

KsM) 5  KA-t) ,  4  > 

where  ac  component  ok)  is  assumed  wi de-sense  stationary  in  the  observation 
interval.  The  latter  term  can  be  developed  as  follows: 

\  ^e[^fr)ey^i2ri4j) 

where  "Rjr)  is  the  correlation  of  the  ac  component  of  the  complex  envelope. 
Now  if 

K  tr)  ~  ^q,  e*r('7lrtr)  05  lr~*  +“' 

h 

each  term  represents  a  pure-tone  component  of  random  phase,  at  frequency 
of  magnitude  Tkl.  If  ^  =  0  f°r  SOme  value  °f  k>  this  term  indicates  a 
pure-tone  component  of  random  phase  at  carrier  frequency  f0. 

Thus  the  deterministic  component  of  the  signal  can  contain  pure- 
tone  components  of  known  phase,  while  the  ac  component  can  contain  pure-tone 
components  of  random  phase.  Coherent  processing  on  both  of  these  components 
is  possible  and  will  be  so  indicated  by  the  optimum  processor. 


(B-19) 


(B-20) 


(B-21 ) 
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APPENDIX  C.  DEFLECTION  CRITERION 


In  this  appendix,  we  consider  maximization  of  the  deflection  at 
the  output  of  the  filter-energy-detector  in  figure  C-l . 


Filter 


F»  |4er  Efi  evTjj  -  D  et  et'tor 


The  input  xlt)  is  composed  of  either  signal-plus-noise  or  noise-alone: 


x/t)- 


sWfhlOl 

or  1  f, r  t**<*»,  T-  V4*, 

nft)  J 


(C-l) 


where  the  signal  and  noise  are  stationary  processes  during  the  gated  interval 
of  time,  independent  of  each  other.  The  only  assumption  on  the  signal  statistics 
is  that  ensemble  average 


Oj  «4beru/>st 


(C-2) 


That  is,  5$  is  stationary  in  the  observation  interval.  A  similar  relationship 
holds  for  the  noise  hffc).  Since  input  x!6  is  a  gated  process,  y ft)  is 


non-stationary.  We  allow  filter  lift)  to  be  unreal izeable,  and  we  Integrate 
the  squarer  output  over  all  time  to  ensure  we  collect  all  the  signal  output. 
The  output  of  figure  C-l  is 
+*•  V*  r  4 


-  £ [j\  x/uhmJ 

=.  J00 cH:  JJ  tlw  <Jv  xfu)x(v)K  H-u)hft-v) 
~  dv  x(w)xfy)  Ah(u-v), 


where  \  is  the  autocorrelation  function  of  the  filter  impulse  response: 

(u-v)  =  Jji  lift  -  u) l»  v) 


Equation  (C-3)  may  be  expressed  alternatively  as 


4  4«*> 

-  fjfiYtf)r  *  j  «iH»rixc 

-Toe  -  °° 


(C-3) 


(C-4) 


(C-5) 


using  Parseval's  theorem,  where 

+r 

H  If )  =  J  dr  e\j>(-i2rrfr)  hk)> 

—  00 

XV) -  £*  e^|>('i2Tr-ft)  x/t).  (c-6) 

Equations  (C-3)  and  (C-5)  indicate  the  generic  form  of  processing  attainable 
by  the  block  diagram  of  figure  C-l,  for  any  time-bandwidth  values  and  any  signal 
and  noise  statistics. 

The  mean  value  of  the  output  of  figure  C-l,  for  signal  and  noise 
present,  is  available  from  (C-3)  and  (C-2)  as 

£^i)5fN]  =  {[  du  dv [Kj (u-v)  +  R„ (u-v)JA  (u-v).  (C-7) 

Then  obviously,  for  noise-alone,  the  mean  output  is 
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£  [?)  w]  =  $  du  dv  K„  (*- v)  Afc  (u"v) . 


(C-8) 


The  difference  in  these  mean  outputs  is  equal  to  the  mean  output  for  signal- 
alone  at  the  input,  and  is 

£^54-Nl}-e^)M}  -  EHs)  =  Jfdudv  T?s(u-vMk(u-v) 


=  { Tck  (t-  hi)??,  Jr)Aj.  fr) 


(C-9) 


where  and  6rs(f)  is  the  signal  power  density  spectrum.  (We  think 

of  underlying  process  5$  as  being  stationary  over  all  time,  but  being  gated 
and  observable  only  during  (taj'HV  If  we  let  be  the  gating  function, 

then  we  have  (more  generally  than  ( C-2 ))  the  input  observable  signal  correlation 

[jft,)  sfe)J  -  all  41;£.)  (C-10) 

Exact  relationship  (C-9)  holds  for  any  interval  T  and  any  signal  spectrum 
(including  a  possible  dc  component,  corresponding  to  a  deterministic  component 
of  sffc))- 

Now  we  make  a  simplifying  approximation:  let 

T»Wj'  Owe)  R5H=0.  (c-ll) 

Thus,  6rs|0  will  have  no  impulse  at  the  origin.  Then  (C-9)  yields  the  useful 
approximation 

-too  4* 

Efrls}  =  Tj  *rfch)\ hr)  =  T  j <Jf  | .  (c-12) 
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Now  we  address  the  noise  output  component.  For  noise-alone  at 
the  input  to  figure  C-l,  we  have 


(C-13) 


Its  mean-square  value  is 


-foO 


(C-14) 


— 


We  now  assume  that  y* (f)  is  a  zero-mean  Gaussian  process;  a  sufficient  condition 
for  this  to  be  true  is  for  noise  input  n(f)  to  be  zero-mean  Gaussian.  Then 
(C-14)  becomes 


-  B  JWfcOSTOjFD  +  2  P UOjq], 

—00 


(C-l 5 ) 


giving 


foe  _  a 

)loir  [a)  n]  =  2  <fe  lit, 


But  consideration  of  figure  C-l  and  (C-l)  gives 


=  if  iw  t)v  (u-v)  Vi  >i  fe  - v) . 


Hence  (C-16),  (C-l 7),  and  (C-4)  yield 


Vor  K)}  s  2  jj  cR,  ft  b  Kw  (u-v)b  fe-u)  life- 


-  z 


vv 

(fj^U  &  \y  Jw  felu-v)fe^-u)AK  (u-/)4  (v-d), 


(C-16) 


(C-l 7 ) 


(C-l 8 ) 
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which  is  exact  for  any  noise  spectrum. 

An  alternative  exact  expression  to  (C-18)  is  available  by  employing 

(C-6) ;  it  is 

MM  =  Jt  ^  |Hlt)f  iHtor  |0Uu)l\ 

where 

5  ft  <Mv  K  (U-^  f- i  2nr  u+  ffc 

is  a  two-dimensional  noise  spectrum. 

For  the  special  case  of  white  noise,  we  have 

R.It>-tMw,  ft, If ) =  -rS"1' 

and  (C-18)  yields 

\  (T 

Var?«|u}  =  iN’JJdwty  Ak(u-;.)  =  _jdr  (t-  W.  (c-22) 

This  expression  is  exact  for  white  noise.  The  alternative  expression  (C-19)  is 
unchanged,  but  (C-20)  specializes  to 

|fl,(f0,-Of  =  (c-23) 

where  5lic(x)  s  5  in  t»T  *)/("■*)• 

If  we  make  the  reasonable  assumption  that  the  frequency  width  of  m 
in  figure  C-l  is  comparable  to  that  of  the  input  signal,  then  for  T1»WS 
and  we  have,  from  (C-22) 

V. rkM  2  T it  KM  -  (C-24) 

'  — —P* 


(C-19) 


(C-20) 


(C-21) 
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the  latter  step  via  (C-4)  and  (C-6).  This  approximation  is  valid  for  large  TV14. 

Finally,  we  define  the  deflection  at  the  output  of  figure  C-l  as 


Vor{?|  N)} 


(C-25) 


Now  calling  on  (C-9),  (C-l 2 ) ,  and  (C-24),  we  obtain  the  deflection  as 

2T 


d* * 


« r^H«r 


(C-26) 


This  expression  is  valid  for  white  noise  and  Its  sensitivity  to  different 

filters  W(f)  in  figure  C-l  is  easily  investigated. 

The  deflection  in  (C-26)  is  maximized  by  the  choice  of  the  Eckart 

filter 


MM  (r  lT  TH. 


The  maximum  value  of  the  deflection  (C-26)  is  then 


d*  =  L  h*  tv%. 

'^0  •— 


Alternative  expressions  to  (C-28)  are 


d° "  i  (hJ  *  TK(w)  ' 


(C-27) 


(C-28) 


(C-29) 


where  the  single-sided  signal  statistical  bandwidth  is  defined  as 


%  = 


(C-30) 
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and  the  received  signal  power  is 

3  =  (C-31 ) 

The  quantity  S/Mc  in  (C-29)  is  the  ratio  of  total  received  signal  power  to  the 
noise  power  in  a  1  Hz  band;  the  latter  quantity  in  (C-29)  is  the  ratio 

of  total  signal  power  to  the  noise  power  in  the  signal  band. 

In  summary,  for  maximum  deflection  at  the  output  of  figure  C-l , 
the  optimum  filter-power-transfer  function  is  proportional  to  the  signal  power 
density  spectrum  (for  white  noise  and  no  deterministic  signal  component);  see 
(C-27).  This  result  obtains  regardless  of  the  size  of  the  input  signal-to-noise 
ratio.  Also  then,  from  (C-4),  (C-6),  and  (C-27), 

Ak>  (t)  =  f  =  k  Ksfr),  (C-32) 

*  _Ho 

which  is  proportional  to  the  signal  autocorrelation  function. 
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APPENDIX  D.  DIFFERENT  FILTEP  CHOICES 
We  have,  from  (45), 


Let  weighting 


C  \ ,  lu)  s  T 
jarVfWy,  |u)>T 


and  define 


|^)\A (v)  R  *  lu  )  6*  p(~  I  2-nrf  U)  *  ^  ^ 


Then  inverting  (D-3),  and  appealing  to  (D-2),  we  have 


(D-2) 


(D-3) 


K .  (u)  -  ]  $  e*p  (i  2irf  u)  (f )  jar  M  K  T>  (D-4 ) 

a  —  m 

Now  we  employ  (D-4)  in  (D-l),  recalling  that  X If)—  D  for 

=  ^»)  1  Jf  (4-^3  ^wlf) 

=  ^jjx  f  df  ^l-f)|xwl  .  (D-5) 


This  relation  holds  true  for  any  weighting  w  satisfying  (D-2). 

Since  the  output  of  the  bottom  channel  of  figure  2  can  be  expressed 


*  i  *rr»r  =  j"  af  iHF)rix(f)r,  w 
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it  would  appear  that  the  choice 


M  • 


(D-7) 


yields  the  desired  quantity  (D-5).  However,  use  of  (D-3)  shows  that 

(D-8) 

where  is  the  Fourier  transform  of  weighting  (D-2).  But  (D-8)  might  be 

negative  for  some  for  example,  if  w(«)  is  chosen  to  be  zero  for  Wl>T 
in  (D-2),  then  V(f)  r  5 in (>rf T)/(ir-f),  which  goes  negative  repeatedly,  and  could 
cause  (D-8)  to  become  negative.  Of  course  the  selection  of  w(u)=l  for  all  u 
yields  VlO1  S  and  (D-8)  is  guaranteed  non-negative;  this  is  the  case 
considered  in  the  main  body  of  the  report. 

Realization  of  Second-Order  Kernel  via  Bush's  Tprhnigma 

Here  we  will  instead  employ  the  technique  given  in  Ref.  7,  pp.  4-6, 
to  derive  the  filtering  block  diagram  necessary  to  realize  X*f-  From  (D-l ) , 
the  second-order  kernel  is  all  with  kernel  transform 

JJ^,  *Ap(-i<0kV  Ra  (V4,)  *  Uw.+  ws)^ak).  (0-9) 


Now  the  most  general  second-order  kernel  is  (Ref.  7,  eq.  7) 


So  if  we  choose 


M  K  <0 


Kc  (fr>)  =  S  fw)  ,  kc  fr)»  )  oil  T, 


(D-10) 


(0-11 ) 


we  have  to  require 


f  (lOr +  (*>*)  k9  (<A>|)  =  £  ^,4  U)z)  Ka  («,)  K>  ('  «,)  *  S  (W'+  ** )K  (*')  (D-l  2 ) 

where  we  have  used  the  real  character  of  impulse  response  )^(r).  That  is. 
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(D-13) 


K.b)KfW  = 

The  block  diagram  of  the  second-order  processor  is  then  indicated  in  figure 
D-l;  some  arbitrariness  is  present  in  (D-13). 


1M.  5  «LCOA«i-  0  ffa  T^te^cr 


For  the  particular  special  case  of  (u>)  *•  K* (<*>),  (D-13)  becomes 

l^<n^)|  =  (D-l  4 ) 

and  we  have  the  block  diagram  presented  in  the  main  body  of  this  report,  when  we 
recall  that  (^(m)  could  contain  an  impulse  at  the  origin  if  ;  see 

(67).  Thus  this  alternative  procedure  yields  the  same  second-order  processor. 


1 
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APPENDIX  E.  GAUSSIAN  FREQUENCY-MODULATION 


The  signal  of  interest  here  is  the  frequency-modulated  process 

sit")  =  A  a»[2irP0l  +  2ir  T  dw ^ (w)  +  <0  j  t K  "h  (e-i  ) 

The  phase  shift  $>is  a  random  variable  and  is  uniformly  distributed  over  2ir. 

The  instantaneous  frequency  of  the  signal  in  (E-l)  is 

■£+fci(0  Ha ,  (E'2) 

where  frequency-modulating  process  *i|t)  is  stationary,  Gaussian,  zero  mean 
(no  loss  of  generality),  and  correlation 

n Ik) *k)  =  W,<  -t,  <E-3) 

The  (double-sided)  power  density  spectrum  of  the  low-pass  frequency-modulating 
process  wffc?  is 


G-*  (f)  =  { c/c  2-n-fr)  Rm  (t), 


( E— 4 ) 


and  the  mean-square  frequency  deviation  is 

B  hi2 If)  =  RJo)  -  (  Jf  Mf).  (E-5) 

The  complex-envelope  signal  is,  from  (E-l), 

Sib)  -  A  2-rj^tlu  wiftf)  +  i  ? 

The  mean  of  this  signal  is  zero,  since  is  uniformly  distributed.  Thus  there 
is  no  deterministic  signal  component. 
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The  autocorrelation  of  *10  is,  for  4«,  *4j  ,4^  <  -tv, 

Rj  ft,  4,)  *  a  Hi) /ft)  =  K  «f[i2^du  ~luj]  >  (E-7) 

where  we  have  assumed  A  non-random.  Since  wi(4)  is  a  Gaussian  random  process, 
the  dimensionless  real  random  variable 

5  2-rrf  da  h,|u)  (E-8) 

is  Gaussian,  with  mean  zero,  and  mean-square  value 

•4*  T  X  x 

y2  -  4 ir  Jf  <lu  K*  (u-v)  *  tir  Jjf  £  <)y  (y)  =  4ir*J  du  (r-  (h) 

-  -  2fj *  IT*  fcX  (E-9) 

where  we  have  employed  (E-3),  (E-4),  and  defined  T  =  |4,-4,).  Therefore  (E-7) 
is  expressible  as 

^5  =  A7  =  k)J,  (E-io) 


using  the  Gaussian  character  of  y  and  (E-9).  Since  the  right-hand  side  of 
(E-10)  is  a  function  only  of  T  =  |4» “4,1,  I  (4)  is  wide-sense  stationary  for 
<.  4  <  4^  ,  and  we  say  that 

Ks  hr)  =  A  z^/c)]  it,  4k  oisevva4  ion  inkrva).  (E-ll) 

According  to  (E-9),  <T*fc)  can  be  interpreted  as  the  variance,  after  T  seconds, 
of  the  "random  walk"  process  y  in  (E-8). 

Since  r7lo)~  0  from  (E-9),  we  find 

\\o)  =  1 5  {Of  "  A2,  (E-12) 
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using  (E-7)  and  (E-ll);  that  is,  the  power  in  the  frequency-modulated  signal  sfO 
is  independent  of  the  frequency  modulating  process  w»|fr)  or  its  spectrum 
Limiting  Behavior 

From  (E-9)  and  (E-5),  we  find 

rT 

(T1  It)  -  ^(r-|u|)^(o)  =  4-rrff^T*,  oil  T>0.  (E-13) 

-T 


Furthermore  (E-9)  yields  directly 


(T  tr)  -  'fir  05  T  -»  0  + 

In  order  to  determine  the  behavior  of  \  hr)  as  we  need 

to  know  the  behavior  of  (r'Kr),  which  in  turn  depends  on  (f)  via  (E-9).  In 
particular,  suppose  that 

(f)  ~  ft  tfl  as  f  0  )  )>  >  -  I 

We  require  i>>-|  so  that  the  mean-square  frequency  deviation  (E-5)  remain 
finite.  Then  we  find,  from  (E-9)  and  (E-15),  the  proportionality  behaviors 


(E-14) 


(E-15) 


X 

<r*  (x)  °<C  ^  l"  x 
c 


•for 

fpr 

V 


)  <  P  <  I 

p  =  |  i  Qs  T~>  +  *>• 
P>  I 


(E-16) 


Thus  for  -l  <  P  <  1?  <T*fc)  tends  to  infinity  as  T  does.  This  in  turn  dictates, 
via  (E-l 1 ) ,  that 


U.M 


Q5  T-»  + ofi- 


0  ,  -\<  P  £  I 

Thus,  according  to  appendix  B,  there  is  no  carrier  term  present  if  -|<«»sl>  but 


(E-l 7) 
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there  Is  a  carrier  term  of  random  phase  if  1)>|,  that  is,  if  faster 

than  -f  at  the  origin.  In  either  case,  there  is  no  carrier  term  of  known  phase 
present,  because  the  uniform  phase  of  $  in  (E-l)  precludes  any  deterministic 
component. 

If  we  assume  that  tr’tr)  approaches  the  constant  value  (T*i*)  as 
then  we  can  express  (E-l  1 )  as 

J?f  (r)  =  (e-is) 

Then  the  spectrum  of  the  complex  envelope  s(t)  is 

£s|f)  =  A*  irlp6)\  5(f)  t  3(f),  ( E-i 9 ) 

where 

Btf)  *  A’  Jcfc  eo<p(r»  M r)[expf-  i  rxlv)}-eKf>[~  i  (E-20 ) 

The  former  component  in  (E-l 9)  represents  a  pure- tone  carrier  term  of  random 
phase  at  frequency  f0|  while  the  latter  term  represents  a  distributed  spectrum 
about  the  carrier;  recall  we  are  dealing  with  complex  envelopes  here.  The 
relations  (E-l 9 J  and  (E-20)  actually  hold  whether  <raW  is  finite  or  infinite; 
in  the  latter  case,  there  is  no  carrier. 

If  6rJo)*0,  then  P  =  0  in  (E-15).  Then  (E-9)  yields 

tMt)  -  4VtJ<*,-RJu)  -  4V^(o)r  «*  T-+«,  ,E-21) 


in  agreement  with  the  first  line  of  (E-l 6 ) .  Thus  in  this  special  case,  there 

is  no  carrier  term,  and  the  signal  correlation  ( E-l 1 )  or  (E-l 8)  decays  exponentially 

with  T- 

Fro*  (£-!),  Here  q)so  fol/OWJ 


(E-22) 


and 


8V 


(E-23) 
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Since  the  integral  on  f“  in  (E-22)  can  be  negative  for  some  values  of  T,  <r*|r) 
is  not  always  monotonically  increasing  with  T,  despite  its  interpretation  below 
(E-9)  as  the  variance  of  the  random  walk  process  y  in  (E-8).  An  example  would 
occur  for  spectrum  )  narrowband  about  a  non-zero  frequency. 

Approximations 

Equations  (E-14)  and  (E-21)  suggest  an  approximation  for  (J*fr) 
when  t^lo)  #0,  namely 


<rl!r) 


4 V  ^  r* 

H- r/fij,  10) 


T2  0. 


(E-24) 


Substitution  in  (E-ll)  then  yields  the  approximation 


"Ks  lr)  A’  exp 


2-rr  <r*  ra 


T>0)  £Jo)£0. 


(E-25) 


We  now  give  several  examples  of  spectrum  GrN(f)  for  the  frequency- 
modulating  process  mitt,  and  illustrate  the  properties  derived  above.  In  these 
examples ,  K  is  a  characteristic  frequency  of  spectrum  M),  but  it  is  not 
necessarily  the  bandwidth  of  Wjftt.  However,  in  the  first  three  examples,  where 
6i,(o)*o.  W.  can  be  interpreted  as  an  equivalent  bandwidth,  in  the  following 
sense: 

=  .h#  ,  AM  =  ,  s-,io)  *  0.  i 

G-*,lo)  Q-„(o)  klo) 


That  is,  K  is  the  width  of  a  rectangle  with  amplitude  equal  to  the  origin  value 
of  the  spectrum  fr-W,  such  that  the  rectangle  and  true  spectrum  have  the  same 
power.  In  this  case,  (E-24)  yields  the  approximation 


4irVwT* 

1-hW^T 


r>0,  #0. 


( E— 27 ) 
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This  approximation  is  plotted  in  figure  E-l  as  the  bottom  curve;  the  remaining 
curves  are  the  exact  results  for  the  first  three  examples  given  below,  and  are 
seen  to  be  well  approximated  by  (E-27). 

Equation  (E-9)  points  out  another  approximation  that  could  be 
extremely  useful;  we  have 

cr'fr)  *  4ttV { $  (r„  (f)  (E~28) 

Now  the  equivalent  width  of  kit)  is  K;  see  (E-26).  On  the  other  hand,  the 
sine*  function  in  (E-28)  reaches  a  null  at  'It.  Then  if 

>  4W*  ,  (E-29) 


(E-28)  yields 

<Tltr)  *  <E-30> 

That  is, 

*  7rr'  -f  t  ■  (E-31) 

But  (E-31)  is  equal  to 

^  -St-  ^  W.-C  -  i  -  IE-32) 

which  is  large  compared  with  unity  if  --“->3.  Therefore  W  Is  substantially 

.  Hi 

zero  for  J>  ,and  we  have 

R,  \r)  =  A’ wr[-  KM)  *  &r  «tl  T,  vf  V  •  ^  >  3.  (E-33) 

That  is,  the  frequency-modulated  signal  correlation  is  approximately  Gaussian  if 
the  RMS  frequency  deviation  is  somewhat  larger  than  the  equivalent  bandwidth 
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of  the  frequency-modulating  process  hjffc).  Some  examples  of  this  comparison  are 
given  below. 

The  spectrum  corresponding  to  (E-33)  Is 


A* 


[-41 


# 


3. 


This  is  a  Gaussian  spectrum  with  standard  deviation  equal  to  <r„,  t  which  Is  the 
RMS  frequency  deviation  of  the  frequency-modulating  process  *nH). 

Examples 
Example  1. 


k)  «  ^  (-2W,  W) 


I  +  (ir-f/W.y 


i  s«,o)  >  on  s.w 

This  curve  Is  plotted  In  figure  E-1 .  ^A)so 


<rak)~ 


05  T-* 


Since  J>  =  0  (see(E-15)  and  ( E-1 6 )X  (J'1h)  =  and  there  is  no  carrier.  The 

distributed  spectral  component  of  the  frequency-modulated  signal  Is  given  by 
(E-20)  as 


3(f)  -  (-» 2 Tpfr) 


The  parameter  D  *  i s  a  measure  of  the  frequency  modulation  Index,  being 
the  ratio  of  the  RMS  frequency  deviation,  Sin,  to  the  equivalent  bandwidth,  WM, 
(characteristic  frequency)  of  the  frequency  modulating  process  *n|fr).  It  is  also 
a  measure  of  the  RMS  phase  deviation  of  the  signal  process,  since  p1  appears  as 
a  factor  in  the  expression  for  <Tl| c)  above. 


-34) 
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Plots  of  are  given  in  Kef.  10,  page  609,  figure  14.4;  however, 

they  are  plotted  in  linear  units  (watts/Hz)  rather  than  in  dB,  and  the  deep 
skirts  of  the  modulated-signal  spectrum  are  not  observable.  We  plot  the  (positive- 
frequency  part  of  the  even)  signal  power  density  spectrum  for  this  example  in 
figure  E-2  in  dB,  where  the  curves  are  all  normalized  to  unit  area  (over 
The  three  parts  of  the  figure  present  the  same  information  on  different  abscissas. 

Figure  E-2(a)  is  most  easily  interpreted  as  considering  W„  fixed,  and  varying 
through  the  frequency-modulation  index  D.  The  curves  labelled  with  an  "E"  are 
the  exact  signal  spectra  obtained  from  the  above  equation  for  B/f),  for  J)  =  2lc, 
Corresponding  to  each  value  of  0  (or  k)  is  also  plotted  the  Gaussian  approximation, 
labelled  with  "G",  afforded  by  saying  that 

<rllr)  *  4V<  t1  ir  all  r; 

although  it  is  recognized  to  be  poor  for  large  x  (see  figure  E-l);  it  can  be  seen 
that  the  Gaussian  spectral  approximation  is  very  poor  in  the  deep  skirt  region  for 
small  values  of  index  D,  but  is  better  for  larger  values  of  index  D.  As  D  increases, 
the  signal  bandwidth  Ws  increases  progressively. 

Figures  E-2(b)  and  E-2(c)  are  most  easily  interpreted  as  considering 
RMS  frequency  deviation  <s^  fixed,  and  varying  modulating  bandwidth  VC  through 
the  index  D.  As  D  increases  (i.e.,  modulating  bandwidth  VC  decreases),  the  deep 
skirts  become  narrower  and  drop  by  3  dB  per  doubling  of  D;  however,  the  skirts 
are  very  slow  in  approaching  the  limiting  Gaussian  spectrum  indicated  by  the 
0  =  «o  curve.  For  J»2,  the  -3  dB  bandwidth  of  the  signal  spectrum  is  virtually 
independent  of  D;  in  fact,  this  is  true  of  any  bandwidth  measure  above  -20  dB. 

However,  it  should  be  noted  from  figures  E-2(b)  and  E-2 (c )  that  the 
-3  dB  bandwidth  decreases  as  D  decreases  below  1  (  VC  i creases  above  «c). 

That  is,  if  RMS  deviation  C  is  held  fixed,  and  the  modulating  bandwidth  Vt,  is 
increased,  the  -3  dB  signal  bandwidth  decreases;  this  result  is  consistent  with 
Ref.  10,  page  609,  figure  14.4b.  It  is  also  consistent  with  the  observation  above 
that  the  RMS  phase  deviation  is  smaller,  since  it  is  proportional  to 
hence  the  frequency-modulated  signal  is  more  nearly  a  carrier  of  small  phase 
deviation  and  a  narrowband  spectrum.  Physically,  the  frequency-modulating  bandwidth 
K  can  get  so  large  that  the  accumulated  phase  perturbation  during  a  time  coherence 
interval  W*  seconds  is  tending  to  zero,  and  the  waveform  is  tending  to  be  progressively 
more  narrowband. 
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W  AW  «is5a  nt>r»n«liid  wi'tli  veif&k  4o  "RMS  •fvep**«cy 

^eVio+i©*  j  S/h«H  D 

Fi^wr^  E"-2.  Niorin«Vi*d  5fan«)  Spctruw;  £*fon«rfo' a |  C&rrtM)  on 
•f^r  Fr*f**»cj  Mo^otimj  JVocwj  wlf). 
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However,  the  figures  clarify  that  the  approach  to  the  narrowband 
limit  is  not  monotonic  at  all  frequencies.  Thus  at  i,  for  example,  the 

spectral  density  initially  increases  as  D  decreases,  until  T>  *  2*  =  \/l6,  at 

which  point  the  spectral  density  decreases  with  further  decreases  in  D.  In 
fact,  by  this  time,  the  normalized  spectral  density  in  figure  E-2(c)  is  well 
approximated  by  the  spectrum 

D 

(f/o-J  +  ’ 

For  f  -0,  this  function  increases  without  limit  as  D  decreases;  for  any  other  ■£ 
it  eventually  decreases  as  D  decreases.  Hence  the  bandwidth  Ws  of  the  signal 
spectrum  tends  to  zero  as  D  decreases,  and  the  assumption  of  large  TWS  will  be 
violated. 

We  will  limit  consideration  to  D>2,  for  which  the  signal  spectrum 
is  approximately  Gaussian,  and  the  bandwidth  on  the  scale  is  virtually 

independent  of  D  (see  figure  E-2 (b ) ) .  That  is,  is  directly  proportional 
to  <TW  for  X)>2.  This  case  also  leads  to  a  significant  RMS  phase  deviation  of 
the  frequency-modulated  signal  in  the  observation  interval  T.  Once  again,  this 
is  in  keeping  with  our  interest  of  detecting  the  most  random  signal. 

Some  analytic  expressions  are  available  for  the  origin  value  and 
the  statistical  bandwidth  of  the  signal  spectrum  3(f).  We  have,  from  above, 

®(o)  '  C  (x  - 1  +  e  ”)] 

where  we  have  used  Ref.  11  ,  3.383  1  and  9.212  1.  For  the  normalized  plots  in 
figure  E-2,  the  quantity  represents  the  origin  value.  A  short  table  follows. 

Notice  that  the  origin  value  increases  without  limit  as  D  approaches  zero,  but 
saturates  as  P->«0- 
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D 

Vp) 

10  lo*  %  (d) 

small 

r 

-10  log  D  -  9.943 

2“6 

6.500 

8.129 

2" 5 

3.273 

5.150 

2‘4 

1.682 

2.259 

2"  3 

.9270 

-0.329 

2'2 

.6035 

-2.193 

2"1 

.4824 

-3.166 

1 

.4364 

-3.601 

2 

.4167 

-3.802 

4 

.4076 

-3.898 

8 

.4032 

-3.945 

16 

.4011 

-3.968 

32 

.4000 

-3.979 

<0 

pF  =  -3989 

-3.991  =  -5  log(fcr) 

Table  E-l .  Origin  Value  of  Spectrum 


The  statistical  bandwidth  of  the  signal  is 

*  „  £jf  Siif)T  .  K;‘io)  r  _j _ 

J(h:?32  hr)  Jdr^/r) 

=  dr  1?  fen*,  M- *  + 1\  p(-2K,lri)]])  ' 

r  2irxJ> _  > 

J»  (l;  2-rr’  p’ i 

by  a  technique  similar  to  that  given  above.  Notice  that  the  statistical  bandwidth 
%  decreases  to  zero  as  0  decreases  (i.e.,  Wm  increases  for  fixed),  as 
noted  above.  Also,  the  bandwidth  on  the  scale  saturates  as  D-*  *•.  A  4*11* 

of  Y4/*w  lions. 
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T> 

HA. 

small 

2-rrP 

1 

3.331 

2-6 

.3070 

2 

3.438 

2-5 

.6  51 J 

6 

3.509 

2-4 

1.148 

10 

3.524 

2“  3 

1.942 

20 

3.534 

2'2 

2.679 

50 

3.541 

2-1 

3.113 

«0 

3.545  *  2i^ 

Table  E-2.  Bandwidth  of  Signal  Spectrum 


A  cautionary  note  is  in  order  here  regarding  the  physical  significance 
of  the  statistical  bandwidth.  Consider  a  unit-power  spectrum 


a,  res+(-^)  +•  ^  £.) 


where 


i  >  M  < 
o,  |u )> 


The  requirement  of  unit  power  means  that  0,4,4-  =  Now  suppose  that 

ox  *  A,  =  I  -  y  o,  -  evpC*'). 


Then 


We  then  find  that 

Q^oO,  A,-*  0,  0,  A,  -*•  Oj 

Qz  -*  \,  I,  0**^  I 


03  V-*  <*• 
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That  Is,  the  spectrum  tends  to  a  flat  one  over  with  power  1;  the 

power  in  the  narrow  component  of  width  a,  tends  to  zero.  However,  the 
statistical  bandwidth  is 


+  L  r 


-x 

e. 


05  X-»  <*>. 


Thus,  this  measure  of  bandwidth  tends  to  zero  even  though  the  spectrum  approaches 
a  flat  one  of  unit  power  over  .  The  situation  is  no  better  for  the  -3  dB 

bandwidth,  for  it  too  approaches  zero  for  this  spectrum,  in  the  limit. 

Calculation  of  Spectrum: 

The  normalized  spectrum  is 


If  we  add 
we  obtain 


6r  (f)  -  [  &*}>(->  2"fc)  ^  hr) 

and  subtract  the  asymptotic  value,  sx^Vb'^-ilj^from  the  integrand, 

r"  =  v-‘;v  +v^  c,s('TrV3W 


where 


3  la)  -  - ')]  «‘ak]'  ’] . 


Using  Simpson's  rule  with  end  correction,  and  sampling  at  intervals  conducive  to 
FFT  processing,  we  obtain,  for  the  integral  above, 

A  J?e  (->  ^h/m)  J  (U) 


where 

jj  (k*)  =  £  (k*)  +  fa  i'Io)  , 
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For  large  0,  this  addition  and  subtraction  procedure  is  not  recommended,  because 
it  requires  the  cancellation  of  two  large  quantities  to  realize  the  result. 
Programs  for  the  spectrum  calculation  are  presented  in  table  E-3. 

Example  2. 

r*)  ,  £„  (f)  -  ’rf'/w;) 

&,(!%)  = . «  SJo) ,  Uk)  =  •  o«  6-Jo), 

J  h)  - 

where 

£(*)  *  • 

The  curve  for  <rft)  is  plotted  in  figure  E-l .  VJe.  <J.*o  l*»ve. 

Jk)  ~  4T*-£jU-  05  T-*  +-• 

Since  P-C,  there  is  no  carrier.  B  (0  would  have  to  be  evaluated  from  (E-20), 
perhaps  via  an  FFT. 

Example  3.  .  OZ? 

KJr)=  l+Kr)1 

Grjl  -  ^  exj>(-2  }fl/K)  ,  t  (^)  OT  &  fo) ,  t  K)  - .  14  fc,  lo) 

if2  Id  =  ft  Qrcfaw^)  "  "a  In  ^ 

where  t/s  The  curve  for  <r*fr)  is  plotted  in  figure  E-l.  A Vso 

<TZhd)  ~  4tt -^-T  05 


There  is  no  carrier. 


3tX3  MSI  is  SCSI  WAXftTT 

PROW  yy'\  fai;  --I;  *,v  jUC 


10  dim  xt  i;£04S'»,Y':i:2048>fc<i:5i3> 

2Q  N=2043  1  M<=2048  N=2'' INTEGER 

30  Pel -.01 
40  PLOTTER  IS  “98721V 
50  i  PLOTTER  IS  “GRAPHICS" 

GO  GRAPHICS 

70  SCALE  6,1 2, -7,0 

80  GRID  1,1 

90  ,  PE HUP 

1G0  S«2--'>-N*4el  > 

11 6  FOR  ,T  =  -i  TO  3 
120  D  =  2-  J 

130  F  =  LGT  (.  SQR  2  +  P  I  ’>*D> 

140  T=. 5- P'2 

150  FOR  1=0  TO  N  2 

1  6  O  U  =  I  *  S 

1  70  V=LGT  EXP  \  —  T  *  U  2  >  )  -F 
180  PLOT  U,Y 

190  IF  V •: - 7  THEN  220 

200  IF  U . 1 2  THEN  220 

210  NEXT  1 

220  PEHUP 

230  NEXT  T 

240  FOR  J=-l  TO  3 

250  B=2  .T 

260  MAT  X*2ER 

270  MAT  Y=2ER 

230  3=1.01*0  1  "'2 

2  9  O  X  (.  1  - .  5 

300  FOR  K= 1  TO  h-1 
310  P=K*De1 

3  2  O  X  <  K  +  1  ■  =  E  X  PC-  3  *  ■  P  1  +  E  X  F'  <  -  P  )  >  ) 

3 3 0  I F  X C! :  +  1  K  1 E - 2 0  THE N  3 5 0 
340  NEXT  K 

FFT  Procedure  in  Viere. 

1050  F=LGT<De 1 > 

1060  8  =  2  N-*  De  1  :■ 

1070  FOR  14=1  TO  N  2+1 
1030  U=S*<:  14-1  > 

1  0  9  0  PLOT  U ,  L  G  T  <  X  i.  1 4  >  +  F 

11O0  IF  U,12  THEN  1120 

1110  NEXT  14 

1120  PENUP 

1130  NEXT  J 

1140  END 


(a)  Program  for  Figure  E-2(a) 


Table  E-3.  Programs  for  Spectrum  Calculation 
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10  GRAPHICS 
20  PLOTTER  IS  "9872A" 

30  V  In  =  6 

40  SCALE  Ct,  Vro,  -6,1 
59  GRII'  .5,1 
SO  PEN UP 
70  DIM  X<  1:2048  >,  Vil  -.2048  >,C  0:51  3  > 

SO  N  =  204  8  !  N<  =  2048  N*2A INTEGER 

SO  COM  S 

100  F= . 5 iLGTi 2*PI > 

110  POR  V=0  TO  Vm  STEP  .05 
l£Ci  V = L G T  < E X P<-.5*V  ' 2 >)-F 
ISO  FLOT  V ,  Y 
.1 40  NEXT  V 
150  PE UUP 

ICO  FOR  Jd-~2  TO  -6  STEP  -1 
170  D=2-  Td 

180  D  e 1  =  . O  2 - D 

iso  mat  >:=zer 

200  MAT  V=2ER 
210  S=  <  p  I  *D  >  -  2 
220  F  -B* 1  D 
230  PoR  r-0  TO  M-l 
240  X<M  1  >*FNG<K*F) 

250  NEXT  Y 

200  output  o : " n  *•*;©;  -Dei  ;x <  l  > ; X <. 2  • ;  x •: 3 > 

2  7  o  0  U  T  p  o'  T  O ;  k  ; : :  <  x  -  2  > ;  x  i  k  —  1  > ;  X  *.  k  1 

2  S  0  X  ■.  1  7  *  X  <  1  T  +  2  ■=  1  *  P  I  2  *  D  *  E  X  P  ‘  S 1  o 

290  0=16  15 

300  FOR  K=2  TO  N  STEP  2 

3  1 0  X Y  '■>  =  X  r-  :>  *  0 
320  NEXT  K 

330  6=14  15 

340  FOR  r  =3  TO  N-l  STEP  2 
350  *  >#Q 

•'CO  NEXT  K 

FFT  Procedure  iw  \ier«. 

10/0  T  =  B*EXF'  ■  S  ' 

1030  D*  1  >•  =  2  'iliHlJ 
1 O90  FOR  V-0  TO  Vm  STEP  .01 
1 100  VI =T  ■ V- 2  +  3  • 

1110  F r- V  !'•=•  1  ’ 

1120  T  = I N I  ■  F r  • 

1  1  i  O  F  r  -  F  r  -  .1 

1140  V 2  =  X  •'  I  +  1  '•  *  •'  1  -  F >  <  +  ,  '.■  .1  +  2  >  * F r 

1  1 5  0  P  L  0  T  V ,  l  G  T  ■.  V  l  ♦  D  s  1  *  V  2  1 

1160  HEX T  V 

1170  PENUP 

1130  NEXT  Id 

! 1 90  END 

1200  DEF  FNu<P> 

1210  COM  S 
1  2  2  0  X  =  S  *  E  X  P  <  -  P  > 

12  30  T  ■=  1 
1240  G=0 

1 750  FOR  H=  1  To  10.00 
1260  T  =  -T ■■  N 
1  270  G  =  G  +  T 

1230  IF  ABS'-T  ■  MBS G  >  *■  1 E- 1  2  THEM  1320 
1290  NEXT  N 

12  00  OUTPUT  0;  "  1000  TERM'S  AT"  IP 
1310  PAUSE 

1370  RETURN  E'  P <  -3*  •.  P-  1  •  •  4-G 
1330  FNEND 

(b)  Program  for  Figure  E-2(c) 

Table  E-3.  Programs  for  Spectrum  Calculation 
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MIS  PAii  IS  BC5I  9UALITT 

PROW  & >y  >  t  j,  *,  ;  .  JH)  ijQfj 


t 

1 


Example  4. 


0  I  v  —3  1-  brYL-r)3  _ 

[i+MUfX 

C„(f)  =  2  «p(-2lf)/vi). 


This  spectrum  goes  to  zero  linearly  at  f=o,  l.e.  )>  *  I  In  (E-15).  It  peaks 
at  we  have 

=  i+  M1*-) ,  t,0'rO  ■  •+'  fr,(*N)> 


r’(r)  *  4  In  [|  +  (ttKtX]  , 


by  use  of  Ref.  11  ,  3.943.  Also 


r'M-sMhr  as  r-»  +  »> 

in  keeping  with  (E-16).  Thus  there  is  no  carrier.  B(f)  follows  from  (E-20)  as 

B(f)=  ~  £[&  c*(2£x)e*r[-24h('W)} 

The  previous  examples  all  had  (Tl^o)  =  «>,  and  hence  no  carrier.  The 
next  two  examples  have  a  carrier,  since  e-*  (Ac*  as  -f-  -*  0  . 

Example  5. 

K,  hr)  =  ()-  Is)  ©xf  (-n-K  T*) 

'  2n-^6-  -^rwp(-T- f‘Ac). 

The  peak  of  6rw(f)  occurs  at  f»  )  then  we  have 
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6U.VC)  =  .31  e»(/.JW.)-.0)4  S-J  5‘H,), 

t*  fc)  -  ^  [l  -  e,<r(-’rWir-)] 

rx  M  *  +ir  -%  +  *■ 

Thus  there  is  a  carrier.  Define  o(  *2Tr^k-  Then 

2-^J^dv  os^i^  x)  -  e -  ^-«r}] 

The  total  signal  spectrum  is  available  from  (E-19)  as 


%(f)=  AV'Stf) +B(f). 


Evaluation  of  B(f)  directly  via  an  FFT  appears  fruitful;  we  have  rapid 
convergence  of  the  integrand  since 


]  ~  I  ~  ^ 


-TTX 


05  X 


Example  6. 


x>  l  \  I  -  3  (irKr) 

K  rr/  -  T  ; - ^ 

'’m  *1 

This  spectrum  peaks  at  -f=WM.  We  find 

tK)=  .s+fcK),  fr,(3H,)=  it 


r1  It)  - 


I  + 


<r*l',)=  4^-  * 
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(T* 

Thus  there  is  a  carrier.  Define  o(  -  2  .  Then 

1  ^  -  **a~  «}] 

*  i  K  e  <i^7  "?T •O*  f0S(X x)  TFMf  • 

Direct  numerical  evaluation  of  Blfi  via  an  FFT  could  take  advantage  of  the 
following: 

<rH  + 


©xi 


\  n  (if?)  * 1  -  -  w-  ■  \T^r'}] , 


and  the  bracketed  term  in  the  last  line  is  OU"*)  tf5 
Validity  of  Approximation  (E-33) 

In  order  to  determine  the  accuracy  of  (E-33),  we  define  the  two 
dimensionless  parameters 

x  =  2-rr<rmT ,  t>  =  • 

Then  we  wish  to  compare  correlations 


*X  ^  [~  a  a*)  toy  [-i**]  , 


for  different  values  of  D.  In  this  normalized  form,  examples  1-3  take  the 
following  respective  forms  for  <r* : 


Wt> *4t±ts - 


wkre  ^  ^  > 

-  ej _ 


J: 


fnr 


-■j  u^ieie  ^ 


fiFS>  * 


£D*[j?x  Qfc-U^  -  £  InG+^x1)]  uk  p 


ip 


(E-35) 

(E-36) 

(E-37A) 

(E-37B) 

(E-37C) 
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Asymptotically,  these  results  take  the  forms 


*’('-53)  05 

2-rrDx-^D1  « 


(E-38A) 


4-J&)  *  x"0+ 


2-rrDx-4TrD  45  X-» 


as  x-*  0+ 


2irDv-8PMo®-?l>1as  x*+-  . 


(E-38B) 


(E-38C) 


Equation  (E-36),  with  each  of  the  examples  in  (E-37)  substituted, 
is  plotted  in  figures  E-3  through  E-5,  along  with  a  good-fitting  Gaussian 
approximation.  We  notice  that  as  D  increases,  the  good-fitting  Gaussian  curve 
approaches  Also  the  approximation  of  (E-33)  is  excellent  for  D>2. 

Thus  when  the  RMS  frequency  deviation  of  the  modulating  process  is  somewhat 
greater  than  the  equivalent  bandwidth  of  the  modulating  process,  the  correlation 
of  the  FM  signal  is  approximately  Gaussian.  See  also  Ref.  10,  Ch.  14. 
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0  .5"  I  1.5  2  2S 

X  *  2r^r 


Figwr*  E~4.  Cor«\«tion5  "Our  £xam^k  2 


APPENDIX  F.  EVALUATION  OF  X a 


We  wish  to  evaluate  as  given  by  (110).  To  this  aim,  we  let 
Un*t*- H*++k)/2»  and  en,P1ov  (113)  t0  exPress»  for  real  ^ 

Tu  -  flff*  r (F1 

Approximate  Evaluation 

Before  we  evaluate  this  quantity,  we  first  approximate  it  by  using 
the  Gaussian  replacement-approximation  mentioned  in  (105)- (.106).  (A  more 
thorough  treatment  is  given  in  (F-42)  et  seq.^T^en  Xu  *5  reflated  by 

Jffi  rk)rk)r^r|i|^  ^ 

O7-2) 

Letting  UK  *  tV+V^  for  the  second  term  in  (F-2)  becomes 

I*  WV’/  *  $(>  (V''  v»V  ^  ^  ^  >  (F-3 


(F— 3) 


where 


*|j«f  rK)r(i4+v,)^kfH)r(^>s).  (F-4) 


For  large  TW„  (F-3)  becomes 


1^  3f  (fv (0)O;  o)  jf[dv,  <Jvx  ^  ^  (vr  w) ^ (vx-V3)^V(p(v,) 


(F-5) 


And  we  have,  from  (F-4)  and  (113),  ^0,0, o)  =  »/V.  Furthermore,  making 
the  Gaussian  assumption  (117)  for  the  signal  correlation  f>,  (F-5)  becomes. 


with  the  aid  of  (125)t 

^JjJ<)v,<Ms  &xp[~ ?W^2  v,V2v’  +3  v£ v,v,-  2v,  v3}]  =  • 


(F-6) 


Employing  this  result  In  (F-2),  we  have 


x  + 


Cntf  ’ 


i.e.,  the  replacement  approximation  leads  to 


(F-7) 


( F— 8 ) 


How  good  this  replacement  is,  quantitatively,  will  be  indicati*  below. 
Exact  Evaluation  for  Large  TW, 

Returning  to  the  exact  expression  (F-l),  let  x  *  u,-i/, ,  y 
v=  U3'i4,  »u  =  i(u3+«t,),  and  then  eliminate  y  via  y  =  V4T;  there  results 


(F-9) 


where  a  shifted  version  of  <j>4  given  in  (F-4).  Now  as  T -*too,  the  ratio 

of  four  p- functions  in  (F-9)  tends  to  1.  Therefore  we  add  and  subtract  1  from 
this  ratio;  the  added  1  leads  to  the  quantity 


jjdx  W{  i  (T) X> v)  '  \j k  ft)]  =  X*  ,  (F-10) 

where  we  employed  (115)  and  (114).  The  remaining  quantity  is 


More  generally,  we  find,  by  use  of  (117)  and  (125),  that  2,  frWs)' '  r ) 

this  includes  (118),  (124),  and  (F-6)  as  special  cases.  When  the  signal  correlation 
is  not  assumed  Gaussian,  depends  on  the  exact  frequency-modulation  statistics; 
an  example  is  presented  at  the  end  of  this  appendix. 
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(F-ll ) 


T  -  T* 
■V  * 


*  If  ^  <fv  i  fr,  y,  *)[#(?,*,*)-)] 

=  4>i.(0)0,o)ffd)r  {?(x)^(v)  J cfc[jf  k,*)  v)-  ij 

-  qjrjj ^  drfytr,x,v)-  0, 


where 


|fc,v)  = 


The  particular  value  of  (F-ll)  depends  on  the  exact  signal 
correlation  considered.  We  will  consider  Gaussian  frequency  modulation 
with  an  exponential  correlation,  as  an  example.  Then  from  appendix  E,  example 
1 ,  we  have 

where  D  -  is  the  ratio  of  RMS  frequency  deviation  to  the  equivalent 

bandwidth  of  the  modulating  process.  For  large  D,  we  have  the  approximation 
to  (F-13) 


2ttV»  r^j 


In  order  to  make  this  limit  agree  with  the  Gaussian  approximation  (117),  we 
set  2v/7^<^=  W5  j  we  then  eliminate  4^,  from  (F-13)  to  obtain 

(Strictly,  this  expression,  in  terms  of  Ws,  holds  only  for  large  D;  however, 
the  dependence  on  D  is  rather  weak,  as  will  be  seen  below,  so  we  use  (F-15) 
for  moderate  D  also. ) 


(F-12) 


(F-13) 


(F-14) 


(F-15) 
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Now  we  can  say  that 


f(vSx)  =  •{-’MS - 1  +  “jHWI  *fM  (M6) 

The  reason  for  this  scale  change  in  (F-16)  is  that  we  can  now  write  the 
convenient  normalization 


(F-17) 


Plots  of  are  given  in  figure  F-l.  It  is  seen  to  be  weakly  dependent  on 

D.  Also,  when  we  make  the  substitutions  *  7=^7  >  V >*»  (p-lli 

we  have  the  clear  dependence  on  TWS: 


(F-18) 


where 


(F-19) 


is  a  function  of  only  the  one  parameter  D,  and 

K(*,v)  -  (F-20) 


We  would  like  to  determine  the  exact  dependence  of  on  D  and 

see  if  it  is  of  the  order  of  1 ,  as  suggested  by  (F-8).  It  should  be  observed 
that  substitution  of  the  limiting  form  (F-17)  directly  in  (F-21)  yields 
<?»(X;*jv)=  ©*j>(xv)  and  hence  a  divergent  integral  on  T  in  (F-20).  Thus, 
the  exact  dependence  of  (F-16)  (or  whatever  example  is  being  considered)  must  be 
preserved  for  accurate  evaluation  o^  (F-19);  this  is  in  contrast  to  the  evaluation 
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of  Jg  and  where  the  Gaussian  approximation  (F-17)  was  applicable;  see 

the  end  of  this  appendix.  The  key  reason  for  the  dichotomy  Is  the  presence 
of  ^-functions  in  the  denominator  of  (F-l),  instead  of  only  In  the  numerator. 

When  a  ^-function  appears  In  the  denominator,  its  exact  rate  of  approach  to 
zero  Is  of  paramount  importance,  whereas  if  it  is  in  the  numerator,  its  rate 
of  approach  to  zero  is  far  less  critical.  Then  when  we  evaluate  any  for 
<^>  0)  it  is  very  example-dependent  and  must  be  handled  very  carefully. 

Numerical  Considerations 

Since  (F-l 9)  must  be  evaluated  numerically,  it  behooves  us  to 
utilize  any  symmetries  involved;  these  symmetries  will  be  derived  for  a  general 
signal  correlation.  First,  since  ^  is  even, 

$>(-*,*,»)  '  *)  ’  Kfov)  -  i] .  (F-22) 


Third, 


Second, 

1o(T>-Y»-v) "  3>h  *>v)  >  K  (-*/  *v/ =  K  (x> v); 

J^v  (v)  K  (x;  v) 

*  £  C*  ft  ft  K(xrv)l 

Third,  (  C  \  \ 

3D(r)>‘rv)  *  ’  K^rV 

^  [$, (w) + 

Fourth,  ^(t,x,v),  >  (<(*)»)  *  K(V)< 

T?  '  W[0^[3pfc».vH  "  2] . 

Now  let  Xsl«,  V=  2y  in  (F-25).  There  follows 

xp  '  [ft  +  »/ 


(F-23) 


(F-24) 


Fourth, 


(F-25) 


- 


(F-26) 
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where 


AM  Pp[t+u-'a) 

(F-26)  is  a  general  result,  applicable  to  any  signal  correlation  example,  where 

foMvfe^r)* 

*  Let  the  result  of  the  integration  on  r  in  (F-26)  be  denoted  by 
and  let  the  result  of  the  integration  on  y  be  denoted  by  3k.  Then 


(F-28) 


where  W}  is  a  set  of  integration  weights  (e.g..  Trapezoidal  or  Simpson), 
ard  Ka  must  be  taken  large  enough  to  include  all  relevant  contributions  to  the 
i ntegral .  Si mi  1 arl y , 

ka  k 

'B(Va)  ^  (2 $A  =  A  (^A)  A  (ka;  JU), 


(F-29) 


and 


MM*)  -  *  '/&-*), 


(F-30) 


where 


Pt>  ((wfk-hl)*) 

A  Y  ((•»-) **A)  ^  (  (»+*-*)  a  ) 


(F-31) 


In  practice,  evaluation  of  (F-30)  for  the  example  of  (F-16)  turned 
out  to  be  exceedingly  slow,  due  to  the  slow  rate  of  decay  of  the  integrand. 
Accordingly,  for  this  particular  example,  a  special  modification  was  employed; 
we  express 


Now  for  Tt  U+y,  the  example  of  (F-l 6 >  yields 

#  -  &*p[-  F  ,  F  -  T>  Jink (^)  J,w ^ Gfj) . 

In  this  case,  then 

4  =  1  F  *f(-  £j]  +  «*p[f  <*f  (-  £)|  “  ^}  • 

u+y  I 

Letting  X*  Fexp(-^)}  and  employing  Ref.  12,  (5.2.4), 

Aa  =  Ti-D^  (e*4  e  *- *)  *  C hi(*)  -In  6- -If]  > 

where 

G-  =  TT  )?  ()-  • 


( F-32 ) 


(F-33) 


(F-34) 


(F-35) 


(F-36) 


By  use  of  Ref.  12,  (5.2.18)  and  (5.1.11),  an  alternative  form  of  A,  is  given 
in  terms  of  the  exponential  integral: 

fy7  =  -n-D [e) tfr) 4  B, (- e)  f  iV  1 2 )» ?r] 

=  4 £,(-£■)  4-2  ^4-2Vj.  (F-37) 

This  latter  form  is  the  one  utilized  here. 
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Hence  (F-32)  becomes 

Afa'j)  «  H.E,(^)+2 

and  (R30)  Is  replaced  by 

A  (yhU)  -  iZjXfafc-* -^(O+fc?  K)*2  U+ It), 


where 


k  =  tV(|-  **?(-£$)■ 


Utilizing  ^  (o>y^  —  0  > q*»«)  A(u;o)~  0>  we  can  combine  (F-28), 
(F-29),  and  0F-39),  to  get  k  . 

4’  yte~2] 

"  "X^E,  (£)+  fe  E,  (-<s)+ 5  )>  %  +  ^), 


There  are  only  two  parameters  to  choose  in  (F-41),  a  and  K.  For  a  given  a, 
we  compute  for  as  large  K  as  necessary  to  realize  an  unchanging  sum;  then  we 
decrease  a  by  a  factor  of  2  and  repeat  the  procedure.  Then  for  the  Trapezoidal 
rule  that  we  used,  since  the  error,  in  approximating  an  integral  by  a  sum, is 
proportional  to  a1,  we  extrapolated  the  approximate  values  to  the  limit  (of  a=o). 
Results  for  this  limiting  procedure,  using  a*.2>.»  ,  .  05"  are  given  in  table  F- 
and  figure  F-2;  the  program  in  BASIC  for  the  HP  9845  follows.  It  is  observed 


that  i 

is  of  the  order 

of  1,  and  increases  linearly  with  D  for  large  1 

JD_ 

To> 

J-t> 

D 

TV> 

J-p 

.25 

.515852 

1 

.418002 

.32 

.416039 

2 

.641552 

.4 

.372120 

3 

.880873 

.5 

.355822 

4 

1 .12359 

.6 

.357359 

5 

1.36758 

.75 

.374109 

6 

1.61220 

Table  F-l . 

Values  of  X©’  for  Exponential 

Correlation 

of  Frequency-Modulating  Process  i*t$ 


(F-38) 

(F-39) 

(F-40) 

(F-41) 
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10 

k=50 

26 

Del-.  1 

30 

It -3 

40 

DIM  R  <■  0 :  408 .) .  F-2  (.1:1 80  > ,  E  < 1 

:  IOC 

50 

REDIN  R  •:  8 :  4*1.  ,R2(l:IO,E(l 

:  K  > 

60 

PRINT  “K  Del  D: “ , K; Be  1 ; D 

70 

COM  T , S 

80 

r = i  cpi*  it  > 

90 

S  =  <  P  I  *  D  )  2 

1  00 

FOR  1' 5 — 0  T 0  4*k 

110 

P'.  Ks  ■  =  F  N  F'  (.  1  s  *  D  e  1  ) 

120 

NEXT  Ks 

120 

T  =  2  *  D  e  1  *  T 

1  4  0 

FOR  Ks-l  TO  k 

150 

R2 1  Ks  —  R 2*Ks  >'-2 

168 

E<Ks>=l-EXP<-Ks*T> 

170 

NEXT  ks 

18c 

81=8 

1  90 

FOR  ls=l  TO  k 

2  0  O 

8  2  =  0 

210 

FOR  L  s  =  1  T  0  k  i- 

228 

T  =  •  R  ■  1  =+Ls>  P'  l  i-Ls.'''  2 

230 

S3= . 5* 1  T+  1  T-2  • 

2  4  0 

T  “R  <  2  *ks  +  2*Ls  '■<  "•.R<2  +  Ks  >*R( 

2  *  L  s 

258 

83  =  33+ . 5*--  T+l  -T  -2> 

260 

T  --  3  *  E  (.  I  s  :■  *E  \  L  =-  > 

2  7  0 

T  =  FHE  1  •:  T :  +FNE  '  ■  -  T ,  +2*  LOG 1 

*  >  ♦ .  i 

260 

33  =  38  -P I  »D*  T.  Del 

2  9  6' 

FOR  Ms*  I  TO  ks+Ls-1 

300 

I  1  =  k  S  +  L  i  - 1'1  i 

3 1 0 

I 2=Ms+Ks  +  Ls 

320 

I  3 =  ft  B  3  (  Ms-  k  s  +  L  s 

3  3  u 

I  4 = M  s  + 1 .  e  -  L  s 

•140 

r = R  •;  11'*  R '  1  2  >  /  ■  R  <  I  3  >  *  R  <  1 4  > 

> 

350 

83=83+ 'T  +  l / T-2  > 

368 

NEXT  Ms 

•3  78 

IF  Ls=Ks  THEN  33=33* .  5 

.  j  6;  O 

S*2  =  82  +  R*; L S  '  * ■Z’  J 

29‘J 

NEXT  Ls 

400 

81=81  +  R  2  i  k  s  1  *  S  2 

4  1  0 

PR I  NT  32  P  I  1 . 5*B*  13*81,  k 

s,R2 

4*  0 

NEXT  ks 

4  30 

STOP 

440 

DEF  FNP < X  ’ 

450 

COM  r,s 

460 

p = ft  £  8  <  x:>*t 

4  70 

r eturn  e ; . p - s *  < p - 1  +  e .■  c p < - p > 

j  > 

480 

FNENIi 

498 

DEF  F ME  1  1  X  >  !  E  .  poner'it  i 

si  It 

Tahle  F-2.  Program  for  (F-41) 


no 


1 


Development  of  Fourth-Order  Correlation 

For  the  signal  process  In  (96),  we  derived  the  fourth-order  correlation 

In  (105): 


i It,) 5%) 5^)2%)  »■  ^’-■9 ^ M)  . 


(F-42) 


If  we  define  the  complex  phase  process 


pit)  =  VX  f  [i  2-r/  cJu  m  |i/)]  ^  4  <  t  <  , 


(F-43) 


then  (96)  and  (77)  enable  us  to  write  (F-42)  as 


where  we  have  made  express  use  of  the  stationarity.  Furthermore,  the  left-hand 
side  of  (F-44)  can  be  written  In  the  two  alternative  forms 


(F-45) 


The  first  form  Indicates  that  If  and  are  kept  fixed,  but  the 

separation  between  these  pairs  of  time  points  is  increased,  then 

the  second  form  indicates  that  if  -t, — t*  and  43-f*  are  kept  fixed,  but' the 

separation  between  these  pairs  of  time  points  is  increased,  then 

Furthermore,  these  are  the  only  ways  these  time  differences  can  increase  without  || 

tending  to  zero.  Thus  we  can  express 

T, (4-  V , V-0  Y4r ^  +  ( 1  + H <  V+’)  >  <F-46> 
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where  the  remainder  function  %  goes  to  zero,  no  matter  how  the  three  difference 
variables  4y-4t;  increase  to  infinity.  Thus  the  first  two  terms 

of  (F-46),  {which  would  be  yielded  for  a  complex  Gaussian  process  fW),  are 
merely  those  terms  which  do  not  necessarily  decay  to  zero  as  the  difference 
variables  increase  in  an  arbitrary  manner. 

Combining  (F-44)  and  (F-46),  we  have  the  very  useful  expression 

which  is  exact,  and  where  tends  to  zero  as  *,1,2  increase,  whatever 

their  fashion.  It  should  be  noted  that  this  expression  and  conclusion  only  hold 
for  a  phase  process  like  that  in  (F-43),  where  ^  is  a  legal  correlation  function 
for  p-  Thus  for  example,  although  |0 W-  ?{-**}*)  is  a  legal  correlation 
function,  it  is  not  a  legal  correlation  function  for  p>  and  R»  need  not  tend 
to  zero  as  its  arguments  increase.  In  fact,  we  find  for  Gaussian 

1U*,  >  (f-48) 

which  does  not  tend  to  zero  for  ^  =  x  + 

If  we  employ  (F-47)  in  (F-l),  there  follows  immediately 

-  £  *  X.  +  Jifa  4|  *  i(x, », 4  />(!)■ -*)%(*, 

s  +  ilpx  Jy  J,, *),  (F-49) 

where  we  employed  (F-4)  et  seq.  Thus  the  first  two  terms  of  (F-49)  corroborate 
(F-2),  (F-7),  and  (F-8);  however,  we  now  see  from  (F-49)  that  the  size  of  the 
remainder  term  is  also  proportional  to  tnwj .  And  indeed,  ( F-l 8)  et  seq.  and 
table  F-l  are  a  numerical  calculation  of  the  sum  of  the  last  two  terms  in  (F-49); 

in  fact,  the  sum  of  the  last  two  terms  in  (F-49)  is  exactly  ( F-l 1 ) : 

>  ft 

Ttx-x*5  (a ~x) (v^ ^ '  (F_50) 


r 1 


Exact  Evaluation  of  Xp  for  Frequency-Modulation  with  Exponential  Correlation 

The  quantities  TV,  Xf.  are  defined  in  (76),  (109),  and  (152) 

respectively.  For  large  "T\^,  these  can  be  simplified  to 

X,  =  Jrj  fc  M 

Xj  -  <*) 

I,  *  r?M  *  *|  (F-! 

If  we  assume  that  signal  correlation  |P  is  Gaussian,  then  we  get  the  approximate 
results  (118),  (124),  and  (F-6).  However,  it  is  of  interest  to  determine  the 
exact  dependence  of  JTj,]  on  0  for  the  example  of  (F-13),  to  determine  the 
accuracy  of  this  approximation.  Making  the  change  of  variable  employed  in  (F-16), 


we  obtain  A 

r. 

- 

D 

(F-! 

T)  =  ^  W> 

r+  -  ^  dw  i>>  fu)^(v-u)^  M^w . 

(F- 

Now  the  question  is:  how  close  are  to  1,  ^4(5^  & 

respectively? 

We  employ  (F-16)  in  (F-53)  to  obtain 

(F- 

By  the  change  of  variable  U*-irP  l*i;there  follows 

ia  =  2-rr*  T>Jo’c)t  i**''  >)). 


(F- 


By  repeated  Integration  by  parts,  this  is  reduced  to 

1  =  (F-56) 

where  ,f;  is  a  confluent  hypergeometric  function.  Evaluation  of  (F-55)  or  (F-56) 
yields  table  F-3.  The  values  of  are  well  tppruimtjaj  iy 


T> 

A 

X, 

.25 

1.32344 

.5 

1.13884 

.75 

1.08798 

1 

1 .06436 

2 

1 .03102 

3 

1.02043 

4 

1.01523 

5 

1.01214 

6 

1 .01010 

7 

1 .00864 

8 

1 .00755 

9 

1.00670 

10 

1.00603 

20 

1 .00300 

30 

1.00200 

100 

1 .00060 

od 

1 .00000 

A 

Table  F-3.  Ta  Dependence  on  D 


1  + 


3-rr^D 


(F-57) 
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which  is  obtainable  from  (F-16)  by  expanding  In  powers  of  1/D.  Thus  a  value 
of  D  =  3  yields  only  a  2 %  error  from  the  Gaussian-approximatlon  value  of  1. 

For  the  evaluation  of  in  (F-53),  we  denote  the  kernel  by 
M  g)V)>  and  note  that 

Kj(-u,-v)=  ki(w>»)> 

Kj(  V,  u)  *  k"j  (U,V),  (F-58) 

to  express  p%> 

-  7^ 

s  ft>  ^ *i ft (l*) ft ft* -ft),  (F-59) 


)J^aWaM 


where  weights  for  k50,  and  z  tor  lfi)»  k>  for  Trapezoidal 

rule,  and  K  must  be  chosen  large  enough  for  negligible  change  in  the  sum.  A 
program  and  numerical  values  for  ( F- 59 )  follow  in  tables  F-4  and  F-5  respectively. 

A 

The  values  of  I3  are  well  approximated  by 


j=_  ,  _  ? 
& 


(F-60) 


A  value  of  D  =  3  yields  a  4.8%  error,  while  a  value  of  D  =  6  yields  a  2.3%  error 
from  the  Gaussian-approximation  value  of 

For  the  evaluation  of  tj.  in  (F-53),  we  denote  the  kernel  by  KM 
and  notice  that 

K+l-ww)  *  kv(w,v,*).  (F-6D 

Therefore 

If  =  ^  ^  (F-62) 
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1 9  B=10 

20  Kc  =  1 00 
30  Eps-lE-6 
40  Del =.2 

50  PRINT  "Kc  Eps  Del  D:  "  ,  Kc  ;  Eps;  De  t 
00  DIM  Pd<. 0:200) 

70  COM  T , S 

80  T=i  <PI*D> 

80  S  =  •:  P  I  * D  >  2 

1 00  F OR  K  =  0  TO  2*Kc 
1  10  P  d  <!  K  >  =  F  N  P  <  K  *  D  e  1  > 

120  HEX T  K 
ISO  $1=0 

140  FOR  K=0  TO  Kc 
l$0  $2=0 

160  FOR  L=-K  TO  t 
170  R=flBS<L> 

)  30  T-Pd  ■:  ft  *  *pd  <  K -L  1 
1$0  IF  h=f  THEN  T  =  T* . 5 
200  $2=$2+T 

210  NEXT  L 

220  IF  K - 0  THEN  $2=32*.$ 

230  Q*Pd<IO*32 

240  S 1=31  +  Q 

250  IF  0  Epi  THEN  280 

260  NEXT  K 

270  PRINT  "END  OF  lCOP" 

280  PRINT  4  FT*Ii«l  24$  1 

290  STOP 

300  IiEF  F NR  i  X  > 

310  COM  T, 3 
320  P=FlB3  (  X  *T 

$30  return  e:;p«:-s*<f-i+e.>:P'  -?>>.> 

340  FNEND 


Table  F-4.  Program  for  (F-59) 


T> 

A 

r. 

T 

1.33412 

2 

1 .23923 

3 

1 .20997 

4 

1 .19576 

5 

1.18736 

6 

1 .18181 

7 

1 .17788 

8 

1.17494 

9 

1.17266 

10 

1.17085 

100 

1  .15630 

QO 

1.15470  « 

Table  F-5. 

A 

Tj  Dependence  on  D 
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using  the  evenness  of 
and  2a  in  v,  we  obtain 


about  X-0. 


If  we  choose  increments  i  in  K 


*4  = 


_!£_ 


)f9{h-V$ 1 

CF-63) 

Li-/.  *  r 

7  j  > 

where  weights  =  £  for  K*0^  and  for  J-k,  and  K  and  L  must  be 

chosen  large  enough  for  negligible  change  in  the  final  sum,A  program  and  numerical 
values  for  (F-63)  follow  in  tables  F-6  and  F-7  respectively.  The  values  of 
are  well  approximated  by 


t,  *  vP 


& 

-JT*P 


(F-64) 


A  value  of  D  =  3  yields  a  7.7 %  error,  while  a  value  of  D  =  6  yields  a  3.8%  error 
from  the  Gaussian-approximation  value  of  V?. 


10  DM  0 
20  K'.  -150 

36  Lr.  *156 
46  E  f:  s  =  1  E  -  9 
56  De i =  .  1 

6  6  P  ft.  INT  "  K  L  c  E  p  &  D  e  i  DS  ",  K  cJlc;E  p  =■  5  D  *  1  1  P 
76  DIM  Pd<..0: 150> 

SO  COM  T.S 
96  r  ~ l .  'PI*  D  > 

1 6  O  3  *  <  P  i  *  D  :>  2 
110  FOR  K  =  0  TO  156 
1  20  Pd'  1-  •  =FUP'  1  *Ii«  i  > 

130  NEXT  t, 

140  S 1 = 0 

150  FOP  1=0  TO  Pc 
160  32=0 

1  7  O  F  0  F  L  -■  I  T  U  L  •: 

1  3  O  T  =  P  d  v  L  *  P  d  H  b  5 ■  L  -  2  *  K  '  .) 

190  IF  L  = K  THEN  T  =  T*.5 
200  52=52+1 

210  IF  T  .  Efji  THEM  230 
226  NEXT  L 
230  52=52  2 

240  IF  K = O  THEN  32=52*. 5 

2  5  0  51=31  +  5 2 

260  PRINT  It  FT  1 . 5*De 1  '  3*5 1 

270  NEXT  1. 

230  5 TOP 

290  DEF  FNP 
3O0  COM  T  ,  5 
310  P=HfcS  <  X  >  * T 

320  RET U R N  E X P  <  - 5 *  <  P- 1  +  E; i P (.  - P  )  )  > 

330  FHEHD 


A 

Table  F-6.  Program  for 
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1 


3 

4 

5 

6 

7 

8 

9 

10 
100 


c6 

Table  F-7.  I* 


r4 

1.78079 
1.58273 
1 .52355 
1.49512 
1.47842 
1.46744 
1.45966 
1.45387 
1 .44938 
1.44581 
1.41733 
1.41421  = 
Dependence  on  D 
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APPENDIX  G.  MEAN  VALUE  OF  'XfiSfN) 

We  wish  to  evaluate  the  mean  value  of  (139)  for  signal  and  noise 
present.  In  this  case,  the  fourth-order  average  of  interest  was  encountered 
in  (94)  and  is  given  by  the  sum  of  (105),  (107),  and  (86).  Substituting  these 
results  in  the  average  of  (139)  yields 

[a+P3M,,VV<-0+2A1  N0fy(v«  S  h-u)  *  sk-th 

A-t)S  +S4,-i)  , 


where 


as  in  (F-44) .  Evaluating  the  delta  function  integrals  and  employing  B-Krh, 


there  results 


"£Vn)  -  -  i  (jz)  I.  -  +  i(£)  (T*  2  r«)  - 


(S-3) 


where 


S  I 


(G-4) 


and  TjjXjj  and  Tfa  are  defined  In  (76),  (109),  and  (110). 

We  have  already  evaluated  1^.  and  ij  in  (118)  and  (124).  We  now  evaluate 
I**- 2Xj,  by  means  of  the  correlation  development  in  (F-47).  First,  from  (F-49), 


Tp-  +  (G-5) 

Next,  employing  (F-44)  and  (F-47),  (G-4)  becomes 

I*.  -  ■^rj/xi' 4, 4 A,  ^  -<)+  tfirlW bV*.f 

=  2iX  +  2i«  +  (6-6) 

where  we  have  also  employed  (76)  and  (F-2),  and  defined 

Q3  =  ^(fS,  +  Y«,-^«r+<)+  <s-7» 

Si  nee  Q>(w)-*0  as  «o  (see  (F-46)  et  seq.)»  we  can  modify  (G-6) 

(by  methods  now  standard)  to 

rf4  =  2  ri  -f  2T^  +  jjff  dx  d*  Q,(v,y,*).  (G-8) 

Combining  (G-5),  (G-7),  and  (G-8),  there  follows 

IK-  2Ttl  =  ^5-jjJ  .  (G_9) 


in  terms  of  remainder  function  defined  by  (F-47).  Since  decays 

to  zero  within  an  interval  of  the  order  of  W*'  in  each  variable,  we  have 


(G-l  0 ) 


where  is  a  quantity  of  the  order  of  unity,  as  in  appendix  F. 

Substituting  (118),  (124),  and  (G-l 0)  in  (G-3),  there  follows 


T*4  ~  t)  (r*J 


-f 


W  frnf 


(G-ll ) 


for  large  TWS. 
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APPENDIX  H.  VARIANCE  OF  T^*(m) 

Let  the  collection  of  correlation  functions  In  (139)  be  denoted 
by  Then  for  noise-alone,  let 

*  j  "•  "I  "j  ",  K«*  >  (H-'  > 

where  an  obvious  shorthand  notation  has  been  adopted.  Then  employing  (86), 

(139),  and  (76),  the  mean  Is 

x"  *  =  t^I(k„u+K,w) 

-  +H1  J  ((■!)■ +  (i-y4  =  -  *  *C T>1*  •  (H-2> 

The  derivation  of  will  be  abbreviated  considerably.  First, 
using  Ref.  6  and  a  technique  similar  to  (H-2), 

X1  r  [  K„.  Km 

=  W.+  J  K  K,jmK„7.,  (h-3> 

I 

where  quantity  A  Is  composed  of  24  delta  function  terms  which  Involve  the  various 
combinations  of  one  even  with  one  odd  number,  out  of  the  numbers  1-8.  For  example, 
one  combination  Is  12,  34,  56,  78;  this  leads  to  the  integral  contribution  In 
(H-3)  containing  KMl3  =  OX- 0  ~  I •  Some  of  the  other  24  terms  are 

far  more  complicated;  for  example,  combination  16,  27,  38,  45  leads  to 

' 2)/Y"  (fofl  "  (H',) 

Substitution  of  all  the  24  terms  In  (H-3),  and  collection  of  like  terms, yl el ds ,  after 
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considerable  labor. 


x*  *  4  M*  +T  ( 1^  -  +  6  *++  3X‘)> 


(H-5) 


where  all  the  terms  Involving  constants,  T^andT,  cancel  out.  Combining 
(H-5)  with  (H-2),  and  recollecting  the  scale  factor  in  (139),  we  get 

Vorfctf’lN  )]  *  i(0  2 


(H-6) 


In  order  to  evaluate  (H-6),  we  appeal  to  the  development  of  the 
fourth-order  correlation  given  in  (F-42)  et  seq.  Specifically,  we  employ  (G-8), 
(F-49),  and  (G-7),  to  obtain 

+61.+ Hi  *  •** ^ " 

[iy  vi,^  -  ^ m eh-*)* y^>(*-3)}  I"-?) 

Alternati vely,  by  use  of  (G-10)  and  (F-50),  we  find 


I?-  2T 

£rws? 


<pi,f 


Finally,  substitution  in  (H-6)  yields 


VmfiTCWi  - 


tptf 


(H-8) 


(H-9) 
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as 


APPENDIX  I.  CORRELATION  OF  Jt^  AND 
The  quantities  Jx  and  are  given  by  (158),  (45),  (77),  and  (139) 

A  -  i  7?  -  Jj*  *  ) 

where  kernel  K  is  the  collection  of  correlation  functions  in  (139)  (as  employed 
in  (140)  and  appendix  H).  For  noise-alone,  the  average  of  the  product  is  (Ref.  6 
and  appendix  H) 

ainXUn)  =  ^  * 

=  3^  J  [\,w  +  fn  Kns  +  K>3»  Y*  +^’s  +/k  ^s7>] 

=  [~ 2^,3  -  + f* fr 3  ~  -  i(^r)  fT»+4rl  (I"3) 

Combining  this  result  with  (71)  and  (140),  there  follows 

i»i>)  -mud  -  -(jjy*,  .  (*-*> 


APPENDIX  J.  EVALUATION  OF 

Is  given  In  (H-8)  In  terms  of  T?p  and  T? ,  which  are  In 


turn  given  by  (F-50)  and  (G-9)-(G-10).  Namely 

I?  *  -  2  it’  + 1  tn (O’  t.  .  (O-D 

Substitution  yields 

tfcf  A,  d*  1?3 +  +(TW,fl+>  (0-2) 

where  two  terms  of  the  -  2  3^  expression  cancel  each  other  by  use  of  the 
property 

if,  (*, »,»)  =  (*.»,*)>  <J-3> 

see  (F-47).  Letting  X-  =  =w/(^K)  In  (J-2),  and 

employing  (F-52),  there  follows 

if  =  JIf  du  clv<W  1?3>  (u,  V,w)  +  * 1*  ,  V-*) 


where 


ILbOw") 5 

A>(v)  A("-»)  l  {  (< 

An  example  of  ^  is  given  in  (F-16). 

Two  useful  properties  follow  immediately  from  (J-5): 

v,w)  - 

They  enable  us  to  express  (J-4)  as 

I?*  ^  ^  ijw  • 


(J-5) 


(J-6) 


(0-7) 
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In  order  to  numerically  evaluate  the  triple-integral  In  (J-7),  It  Is  nec..  'ary 
to  determine  the  regions  of  significant  value  of  ft^,and  to  terminate  the 
Integrals.  Detailed  evaluation  of  7?33  for  the  example  of  (F-16)  for  D  ■  3 
was  undertaken,  with  the  result  that  the  first  term  In  (J-7)  can  be  well-approxi¬ 
mated  by 


where 

L(v)  *  6.5 -  .  0*V,  V--  +0. 


(J-9) 


Letting  the  integral  on  u.  be  denoted  as  and  the  Integral  on  V  denoted 

as  B(»>,  we  have 


Also 


and 


«•»«£»,  m 

A  (ty")  =  J  ^  V*  (u'v>w)  "  ^>P  (-l  ^ » Y>  W), 

— lK) 

M  a,  *  M,’h[w,Uv)}  ♦  L(v);  •  l  f'  . 


(J-ll) 


(J-12) 


A  program  for  the  evaluation  of  (J-9)-(J-12)  Is  given  in  table  J-l .  Extreme  care 
Is  necessary  In  the  evaluation  of  fljp,  due  to  the  presence  of  ^  functions  In 
the  denominator  of  (J-5).  Exercise  of  the  program  for  sampling  Increments  .5 
and  .25  yielded  the  value  Q  *  5.75415  for  D  =  3. 

Combining  this  result  with  table  F-7,  we  find 

II.  f+«  -fir  t>  =  3.  <J-'3> 
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It'  Vc=50 

20  D  e  I  = .  2 5 

30  D«3 

4  0  K  1  =  1  c  P I  *  D  > 

r.u  F2='..Pi^ro  2 

CO  COM  Lei , L ' , F 1 , F2 
70  Kc  -  1 NT  t  Vc  ■'  Del  *41 

SO  Dw*Vc -Kc 
90  9=0 

100  FOR  K-0  TO  Kc 
lie  r=Fni.-:K 
170  IK  K — O  THEN  T=T * . 5 
130  S=£+T 

140  PRINT  4 P I  ’  1 . 5 4 D m * S .  K 

150  NEXT  K 

ICO  END 

176  DEF  F  N P 1 V  1 

ISO  COM  Del  ,  Lv> 

190  Lv**6.5-.08*V 

2‘  0  O  L_ '_  =  I  N  T  ■  2  ■*  L ■  Del  '•  + 1 

2 1  O  D  w = 2  *■  L  v  L  c 

220  W0  =  V-I_m 
230  S=0 

240  FOR  L=0  TO  Lc 

c:  5  0  S  ®  S + F  N  ft  V ,  W  0  +  L  -*■  D  w  > 

260  NEXT  L 

2  i  ‘  0  RET  U  R  N  P  m  *■  S 
230  FNENP 

290  DEF  FNrt  <:  V  ,  U  > 

300  COM  Del, Lm 

310  U0=HIN<  W,  L w 

3 2 0  M .:»I  NT  <  <  UO  +  L >  / D e  1  ■  -*- 1 

3  20  D '.<=•:  OO+L'  O  - Me 
3  4  0  S  =  O 

350  FOR  M ■=  O  TO  Me 
3  6  0  T  -  F  N  R  2  <  -  L  c  +  M  %  D  u ,  V ,  U  1  2 
370  IF  M-Mc  THEN  T=T*.5 

3  3  O  S  *  '3  +  T 
390  NEXT  M 
40O  RETURN  Du 
410  FNENP 

420  PEF  FNROm.I,  V,W  > 

4  3  0  H  1  a  F  N  G  IJ  ■  +  F  N  G  v  W 
440  ft  2  =  F  N  G  •  W  .'  +  F  N  G  ■  V  -  U  ■ 

450  ft?  =  FNG <  V  ■  +FNG <  W-IJ  > 

460  RETURN  Ei  RC-ftl-ft2  +  ft3 '-EXPC-ftl  '-F 

470  FNENP 

480  DEF  FNGOO 

4  9  ( '  COM  P  e 1 , L " , E 1 , F  2 

500  P«flBS<X>*Fl 

510  RETURN  F2*  t  P -  1 +EXP < -P >  > 

520  FNENP 


Table  J-l.  Program  for  (J-9)-(J-12) 
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In  equations  (153),  (154),  (156),  (157),  and  (162),  the  quantity 

£  2.U2  £r  »*3) 

appears;  It  Is  of  the  order  of  unity,  as  anticipated. 

Finally,  reference  to  (J-1),  (F-52),  table  F-l ,  and  table  F-7 
enables  us  to  evaluate 

r*  =  4. 4^  ir  J)  *  3. 

In  equations  (146),  (148),  (149),  (156),  (157),  and  (159),  the  quantity 

(=  i.m  £>r  J>  =  3) 

appears;  it  is  of  the  order  of  unity,  as  anticipated. 


(J-1 4) 

(J-1 5 ) 

(a-  \e) 
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APPENDIX  K.  DERIVATION  OF  AND  ITS  MEAN 

The  cumulants  through  order  4  were  given  in  (39).  The  fifth  (and 
all  odd-orders)  cumulant  is  zero,  as  shown  in  (1 31 )- (1 33 ) .  The  sixth  cumulant 
is  given  by  (Ref.  4,  3.43) 

--  <r‘>  -is<Y+><V>  +  30<>t>  .  <k->> 

From  (131 )  and  (132), 

<rV  #-M‘  -1T|CI‘.  <K'2) 

where 

C  -=^-jxp’  (K-3) 

Here  pl$  is  the  complex  phasor  process 

ylfc)  =  evj[l“2rJ^  cJum(w)]  }  ^  ^  .  <K-4) 

Using  relations  similar  to  (K-2)  for  lower-orders,  we  find  from 

(K-l) 

XT=  jf  -  1  <|CI+><)C|“>  +  12  <10*/)  ;  (*-S) 

where  the  averages  are  over  the  frequency- modulating  process  toft)  in  (K-4), 
or  equivalently  over  the  phasor  process  pit).  Now  from  (K-3), 

<'cO  = 

<  ici‘> = (0  H  *•  •• y*  kMS(  <*-«> 

Ktt=<3>rr.> 

K,jm=  >  Kmm  " 


where 


Substitution  of  (K-6)  In  (K-5)  yields  the  desired  result 

7^  9  \K  a**  “  ^  Ki  Xy*Kf\  .  ( K-8 ) 

Evaluation  of  the  kernel  of  (K-8)  would  require  evaluation  of  the 
sixth-order  average  in  (K-7);  although  possible,  this  is  very  tedious  and  sheds 
no  light  on  the  processor.  Also,  we  can  evaluate  the  mean  value  of 
without  knowing  all  the  details  of  the  kernel  of  (K-8).  Using  the  properties  of 
complex  white  Gaussian  noise,  we  have 

-  if  (£)  NTKJ  Otw+kw*  Ksj»s>  Kn>*» 


-  ^  +  K)l3y  k)W(  iCjy-h  k)S-3)  Kw  Kn +  k,y»  K>] 

•+■  )2  ^  K„  k5s  +  k,  Vjy  K55  +  Kq  Vj,  Y>ss  ^  KsKai  ksa  ^ 

+  KK»Ks,  +  KXX,}\ 

+  )2{)-»- 3^+213]] 

■  +80feJ  Tj  ■ 

Tliertfore 

^  =  TfeM  *  ^  *r  W  ™- 


(K-9) 


(K-10) 
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